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ABSTRACT 
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Hhen  designing  an  active  sonar  homing  torpedo,  certain 
operational  torpedo  parameters  such  as  speed,  turn 
rate,  etc.  have  to  be  decided  upon.  For  a given 
homing  torpedo,  there  must  exist  tactical  guidelines 
of  how  to  employ  the  torpedo,  i.e.  which  firing 
position  gives  the  best  chance  of  a hit.  This  thesis 
attempts  to  gain  some  insight  into  the  detection 
process  during  the  torpedo  run,  as  well  as  getting 
some  indications  of  the  relative  importance  of  the 
different  torpedo  parameters  and  the  tactical 
situations.  A simulation  model  was  used  in  order  to 
generate  the  data  base  for  analysis.  The  results 
stress  the  importance  of  a good  firing  position  as 
well  as  show  how  it  is  possible  to  counter  a bad 
firing  position  by  a high  speed  torpedo.  They  also 
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The  following  analysis  examines  the  performance  of  a 
homing  torpedo  against  a surface  ship.  A homing  torpedo  is 
described  as  a torpedo  which  is  searching/snaking  on  each 
side  of  its  main  course.  It  is  searching  for  a target  by 
transmitting  with  its  sonar  and  listning  for  an  echo.  Ref. 
Fig.  1.  Passive  searching  torpedoes  and  homing  torpedoes 
going  in  circles  are  not  investigated  in  this  paper. 

The  torpedo's  performance  is  a function  of  many 
variables.  These  variables  are  divid  .d  into  two  groups; 
technical  variables;  speed,  max  torpedo  run,  sweep 
angle,  technical  detection  range,  lobe  characteristics 
and  turn  rate. 

tactical  variables  ; firing  range,  attack  angle, 
target  speed,  type  of  target  and  tactical  detection 
range  (sonar  conditions) . 

No  attempt  is  made  to  analyze  the  first  grcup  of 
variables;  instead,  technical  variables  used  are  those  of 
present  technology.  We  are  assuming  a 'standard  homing 
torpedo'  based  upon  homing  torpedoes  in  operational  use 
today  [8].  This  'standard  torpedo'  assumes  conventional 
warhead  and  active  sonar  transmission  for  detection.,  and  is 
unguided . 
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The  technical  variables  (torpedo  parameters)  are  ir.  many 
ways  interrelated.  Por  example,  the  maximum  detection  range 
will  determine  the  transmission  rate,  since  the  transmitted 
energy  must  have  time  to  traverse  out  to  maximum  detection 
range  and  return  as  an  echo  before  the  next  transmission,  at 
least  during  tne  search  phase. 

it  the  same  time  the  torpedo  is  transmitting,  it  is 
searching  (changing  course)  for  a target.  In  each 
transmission,  the  transmitted  energy  is  focused  within  a 
narrow  beam  (lobe) . During  reception,  the  echo  is  confined 
within  the  same  narrow  beam(lobe).  Concurrently,  in  the 
time  between  two  transmissions  the  turning  rate  of  the 
torpedo  must  be  limited  to  ensure  that  the  receiving  lcbe  is 
not  outside  the  direction  from  where  an  echo  may  return. 
Thus  turn  rate  should  be  a function  of  detection  range  and 
the  lobe  pattern. 

In  order  to  maintain  torpedo  speed,  the  number  of 
degrees  of  sweep  on  each  side  of  the  main  course  must  be 
small.  If  the  sweepangle  is  small,  however,  the  width  of 
the  possible  detection  lane  will  be  small  as  well,  and 
consequently  the  detection  probability  might  be  reduced 
during  transit.  Also,  a high  torpedo  speed  creates  a great 
change  in  torpedo  position  between  each  transmission.  In 
this  way  the  torpedo  may  scan  outside  a target  in  the  sweep 
lane.  In  other  words,  the  coverage  density  of  the  lobe  may 
be  low  as  a result  of  the  high  movement  rate. 

As  we  recognize  the  relationship  between  torpedo 
parameters,  tactical  variables  and  torpedo  performance,  we 
know  that  frequently  within  the  naval  establishment 
decisions  have  to  be  made  with  regard  to  torpedo  parameters 
and  tactical  doctrines.  In  localizing  and  defining  these 
relationships  this  analysis  may  be  a tool  in  this 
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decisionprocess. 


The  measure  of  effectiveness  by  which  different 
alternatives  will  be  judged  will  be  detection  probability, 
by  which  is  meant  the  probability  that  the  torpedo's  active 
sonar  detects  and  begins  to  track  the  target.  This 
probability  will  be  measured  by  a digital  computer 
simulation,  construction  of  which  was  a major  part  of  the 
author's  effort  in  writing  this  thesis. 
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II.  NATURE  OF  THE  PROBLEM 


A.  DEFINITIONS 


Lobe  width  is  the  number  of  degrees  fron  the 
centerheading  of  the  torpedo,  until  the  first  minimum  in 
transmission  intensity  is  reached.  See  Fig.  4. 

Detection  range  is  the  range  .to  the  target  when 
detection  first  occurs. 

Technical  detection  range  is  the  max  detection  range 
which  is  technically  and  reasonably  possible  considering 
power  transmitted  and  lobewidth.  It  is  the  basis  for 
determining  the  transmission  rate. 

Aspect  is  the  angle  measured  from  the  positive  direction 
of  the  longitudinal  axis  of  the  target  to  a line  joining  the 
centers  of  gravity  of  the  target  and  the  torpedo. 

Attack  angle  is  the  aspect  at  the  start  of  the  torpedo 

run. 


Sweep  angle  is  the  maximum  number  of  degrees  the  torpedo 
will  turn  off  the  main  course  during  search. 


All  dimensions  are  in  meter,  second,  meter  per  second, 
degree,  degree  per  second.  Speed  of  the  target  and  the 
torpedo  are,  however,  always  given  in  knots. 
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It  is  assumed  that  all  firings  are  successful,  and  the 
torpedo  Mill  not  deviate  from  its  ordered/calculated  course 
and  speed. 

All  firings  are  made  with  a deflection  angle;  i.e.  the 
torpedo  is  given  a course  to  a predicted  hitting  point  with 
the  target. 


B.  ASSUMPTIONS 


Not  only  in  order  to  keep  the  problem  tractable,  but 
also  because  of  modern  torpedo  development,  only  surface 
targets  are  considered.  Previously  within  NATO,  torpedo 
developments  seemed  to  start  as  a development  of  an 
anti-submarine  torpedo  with  later  modifications  in  order  to 
make  the  torpedo  dual  purpose.  However,  today  there  are 
some  indications  that  the  anti-surface  ship  requirement  is 
coming  into  the  development  early  in  the  planning  process 
[7;8;9].  The  entire  problem  is  then  kept  in  two  dimensions. 
The  vertical  axis  is  not  significant  as  we  assume 
isovelocity  condition,  and  we  assume  for  simplicity  that  we 
have  negligible  surface  effect. 

Also,  if  the  intensity  of  the  echo  is  above  detection 
threshold  level,  the  target  is  detected  with  probability 
one.  Probability  of  false  contact  is  assumed  to  be  zero. 

The  main  purpose  of  a homing  torpedo  is  to  counter 
uncertainty  in  target  data  at  firing  and  target  maneuvering 
after  firing. 

For  simplicity  the  following  assumptions  are  made: 

the  target  remains  on  a steady  course  after  firing, 
estimated  target  data  is  used  in  solving  the 
deflection  angle  problem 
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TAM  = target  estimated  speed 

ASP  = estimated  aspect 


= (target  estimated  coarse)  - 

(bearing  to  torpedo) 

TO  = torpedo  speed. 

See  Fig.  2. 

The  difference  between  the  target  data  and  the  target 
estimated  data  are  defined  as  errors  in  the  target  data. 
These  errors  are  assumed  to  be  random  variates  and  are  given 
as ; 

target  range  error  is  uniformly  distributed  between 

- 15  % and  ♦ 15  % of  actual  target  range 

target  course  error  is  uniformly  distributed  between 

- 15  and  ♦ 15  degrees 

target  speed  error  is  normally  distributed  with  mean 

0 and  standard  deviation  3 knots. 

These  errors  are  assumed  to  cover  errors  in  the  fire 
control  solution  at  the  time  of  firing  as  well  as 
non-radical  maneuvering  of  the  target  during  the  torpedo 
run. 

As  shown  in  Eg.  2.1,  estimated  range  does  not  enter  into 
the  calculation.  Estimated  range  would  only  be  used  for 
some  more  complicated  tactical  situations  as  angled  torpedo 
firing  off  the  firing  coarse  of  the  firing  unit.  These 
situations  are  not  covered  in  this  study. 
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III.  PROBLEM  SOLVING  APPROACH 


At  firing,  the  initial  course  of  the  torpedo  is 
uniformly  distributed  between  minimum  course  and  maximum 
course  (main  course  ♦/-  a fraction  of  sweep  angle)  . 

Immediately  after  firing,  the  torpedo  starts  'snaking*. 
During  snaking,  the  torpedo  is  continously  changing  course 
left  or  right  out  to  the  given  sweepangle,  then  back  past 
main  course  and  out  to  sweepangle  on  the  other  side  and  so 
on.  The  torpedo  is  turning  with  the  given  turnrate.  Luring 
the  whole  process,  the  torpedo  is  also  transmitting  and 
listening.  Transmission  interval (TTIME)  is  given  by 
technical  detection  range  as; 

TTIME  = 2 x TEDEC/1 500  seconds  (3.1) 

where 

TEDEC  = technical  detection  range  in  meters. 

1500  = speed  of  sound  in  salt  water,  m/sec. 

The  torpedo  run  is  conducted  in  steps.  Every  0.5 
seconds  interval,  all  positions  and  courses  are  updated. 

At  each  transmission;  the  relative  bearing  to  target, 
and  the  target  aspect  are  calculated  in  order  to  establish 
the  intensity  of  the  echo. 

When  a detection  occurs,  the  following  data  is  stored; 
detection  range  to  the  center  of  the  target. 
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detection  range  to  the  nearest  part  of  the  target, 
detection  bearing  (relative)  to  the  center  of  the 
target. 

detection  bearing (relative)  to  the  nearest  part  cf 
the  target. 

- target  aspect. 

In  addition  to  the  detection  probability,  the  range  at 
which  the  detection  first  occurs  is  also  of  interest. 
Therefore,  we  store  these  data  at  the  first  detection. 

However,  successive  detections  are  also  important.  As 
part  of  the  criterion  for  the  decision  of  when  to  go  from 
search-phase  to  attack-phase,  the  number  of  successive 
detections  (with  no  non-detection  between)  may  be  employed. 
In  real  life  there  is  always  a positive  probability  of  false 
detection.  Sven  if  we  are  not  addressing  the  problem  of 
false  contact  as  such,  we  can  cover  the  possibility  by 
requiring  the  torpedo  to  have  at  least  two  successive 
detections  before  going  into  attack-phase.  Accordingly,  we 
store  also  the  previously  listed  data  at  the  second 
successive  detect  ion (two  immediately  following  detections), 
at  the  third  and  so  on,  up  to  and  including  5 successive 
detections.  This  listing  of  detections  will  give  an 
indication  of  the  decrease  in  detection  probability  if  a 
large  number  of  successive  detections  before  going  into 
attack-phase  is  required  in  order  to  decrease  the 
probability  of  false  contacts. 
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IV. 


MODEL 


A.  SEARCH 


For  simulating  the  torpedo  search,  a Fortran  IV 
simulation  program  was  developed. 

The  program  was  divided  into; 

Main  program,  including  generation  of  statistics  and 
print  out  of  summary  after  all  the  runs  were  completed. 
Subroutine  PARMET  for  setting  tactical  situation  and 
torpedo  parameters. 

Subroutine  FIRING  which  calculates  estimated  target 
data,  and  the  deflection  angle. 

Subroutine  POSIS  which  calculates  the  torpedo  course, 
and  torpedo  and  target  positions  at  each  time  step. 

- Subroutine  DETECT  which  checks  if  the  target  is 
detected  and  if  so,  store  detection  data. 

See  Fig.  3. 

See  Appendix  A and  Appendix  B. 
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B.  DETECTION  MODEL 


A contact  occurs  when  the  acoustic  energy-pulse 
generated  at  the  transducer  and  reflected  from  the  target  as 
an  echo,  is  at  or  above  threshold  level.  In  the  following 
discussion  we  assume  that  the  contact  meets  the  tactical 
requirement,  and  accordingly  we  use  the  term  detection. 


1 • Pet ection  Threshold 


Deciding  if  a detection  occurs  is  a function  of 
detection  threshold (signal  to  noise  ratio),  the  range  to  the 
target,  the  target  strength  and  the  relative  bearing  to  the 
target,  given  a level  of  radiated  intensity. 

The  detection  threshold  for  a torpedo  is  a function 
of  design  and  technological  sophistication  of  the  torpedo. 
Without  making  any  assumption  about  these  variables  in  the 
model,  we  start  with  a given  technical  detection  range,  a 
'standard'  target,  and  calculate  intensity  of  echo  at  that 
range  for  target  aspect  equal  to  90  degrees  (maximum  target 
strength)  and  relative  bearing  to  the  target  equal  to  zero 
degrees.  This  echo  intensity  is  then  the  detection 
threshold  for  every  transmission  during  a run.  If  any  echo 
intensity  is  above  the  detection  threshold,  it  is  detected; 
if  below  the  detection  threshold,  it  is  not  detected. 


2 . Echo  Intensity 


In  calculating  echo  intensity  we  must  separately 
investigate  the  important  factors,  which  are  transducer 
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a.  Lobe 

Characteristics 

The 

transducer  has  a 

main  lobe  and 

many 

sidelobes  as  a 

function  of  the 

transducer's  gain 

and 

relative  bearing. 

Uric*  [ 6 ; 51-57] 

discusses  some  of 

the 

different  types  of  beam  pattern  (lobes) , and  the  following 
mathematical  model  was  developed  and  found  to  give  an 
acceptable  pattern; 


G(9) 


G 

0 


SIN  (x  IT)  | 

1 COS  (6/2) 


X TT 


(4.1) 


where 

X = 9/0 

0 

and 

G = maximum  gain 
0 

0 = relative  bearing 

6 = lobe  width. 

0 


(4.2) 


This  model  will  produce  the  gain-pattern  as 
shown  in  fig.  4. 
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Figure  4 - DISTRIBUTION  OF  LOBES  AND  INTENSITY 
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b.  Reduction  in  Intensity  due  to  Range 

Primarly,  the  reduction  is  due  to  two  effects; 
spherical  spreading  and  absorption. 

Spherical  spreading  is  a known  function,  but  absorption  is 
dependent  upon  transmission  frequences,  water,  salinity  etc. 
In  order  to  simplify  the  model  and  since  spherical  spreading 
has  the  greatest  effect,  only  the  spherical  spreading  for 
reduction  in  intensity  is  considered.  This  reduction  in  a 
one  way  propagation  is  given  by; 


where 


I = I /R 
0 


= radiated  intensity  at  one  meter 


(h-3) 


I = intensity  at  range  R. 

R = range  in  meters. 

c.  Target  Strength  and  Target  Aspect 

When  the  transmitted  energy  pulse  hits  the 
target,  some  of  the  energy  is  reflected  back  to  the 
transducer.  The  echo  intensity  is  a function  of  the  shape 
and  dimension  of  the  target,  type  of  reflective  material  and 
aspect. 

It  should  be  noted  that  the  notion  of  target 
strength  represents  the  ratio  between  target  cross  section 
and  the  surface  of  a sphere  of  radius  1 meter,  or  if  in  dE, 
10  times  the  log  of  this  ratio;  base  10.  In  most 
references,  the  target  strength  or  the  target  cross  section 
is  given  abeam  of  the  target,  see  [ 5 ; 97  ], [ 6 ; 274  ],  without 
presenting  the  cross  section  as  a function  of  the  aspect. 
Urick  [6;  282,283]  gives  , however,  as  figures,  an  indication 
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of  how  the  target  strength  (in  dB)  varies  with  the  aspect. 
Cox[3;60]  states  that  it  will  vary  between  10  and  25  dB. 
All  Measurements  in  dB  in  the  two  references  are  relative  to 
1 yard  as  unit  for  range.  0rick[ 6 ; 283-286  ] indicates  that 
his  reference  (as  given  in  the  figures)  will  not  change  in 
any  considerable  degree  with  changes  in  frequencies  (20-60 
KHz)  or  for  different  targets  (submarines/sur face  ships). 

Assuming  a torpedo  with  transmitting  frequency 
between  50  and  60  KHz,  we  get  a wavelength  varying  between 
2.5  and  3.0  ca(0.025  - 0.03  meters).  As  any  reflection  from 
a target  is  mostly  determined  by  target  fora,  size,  aspect 
and  wavelength,  we  may  use  a model  from  radar  theory  in  our 
next  step.  The  justification  for  this  use  is  that  in  radar 
theory  we  are  wording  in  the  same  area  of  wavelength  and 
target  dimension  as  an  active  sonar  for  a homing  torpedo. 

Crispin  and  Siegel  [4;86]  give  for  target  cross 
section  a model  for  an  ellipsoid  where  the  incident 
angle (target  aspect)  is  a variable.  The  relationship  is  as 
follows; 


6 


2 

IT  • a 


2 


•c 


2 2 2 2 ' 2 Z 2 2~~  2 

(a  -SIN  0-COS  f>  ♦ b -SIN  0 -SIN  9 ♦ c -COS  0) 

(4.4) 


a,  b,  c being  half  axes  of  the  ellipsoid. 
3ef.  Pig.  5. 
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As  we  assume  that  the  transmitting  pulse  is 
always  in  the  horizontal  plane,  9 is  90  degrees,  which  gives 


2 2 2 
lT*a  -b  *c 

— 2 2 2 2 7~ 

(a  -COS  1>  ♦ b -SIN  1)) 


(*♦.5) 


<|  = aspect. 

Urick  [6;275]  gives  for  target  section  a model 
for  abeam  or  ahead  cases,  which  is; 

2 

t = <3/  (41T)  = (b  -c/2  - a)  , 

identical  with  Eg.  4.5.  Note  that  Eg.  4.5  is  an 

expression  for  the  target  cross  section. 

Haslett  [5; 139]  gives  for  the  target  cross  section  a model 
for  both  ahead  and  abeam  cases.  His  model  eguals  Eg.  4.5 

2 

times  a factor  R , where  R is  acoustic  reflectivity 
coefficient  (per  cent)  = 94. 

The  advantage  of  using  Eg.  4.5  is  that  it  gives  the  target 
cross  area  as  a continous  function  of  the  target  aspect. 
Por  our  model  we  will  only  use  the  lower  part  of  the 

ellipsoid  to  simulate  the  ship  hull  below  the  water  line. 

Combining  Eg.  4.5  and  acoustic  reflectivity 
coefficient  we  get  the  following  model  for  the  target  cross 
section ; 


2 2 2 2 
iT-a  -b  -c  • R 

— 2 2 2 2 2~ 

(a  -COS  0 ♦ b -SIN  p) 


(4.6) 


With  reference  to  Urick's  figures  [6;283]  where 
the  pattern  of  the  target  strength  is  given  as  a function  of 


W 


aspect  in  Pig.  9.13,  and  reproduced  in  this  analysis  as 
Pig.  6.  a,  we  still  have  not  obtained  a iiodel  which  gives 
the  same  type  of  pattern.  By  applying  the  following  scaling 
factor  to  Eg.  4.6  we  have  approximated  his  information: 

2 

0 = (0.251  635  •()  - 0.18555  • <p  ♦ 0 . 036  5 • SIN  (3  • ((>  ♦ 0.  17453)) 

2 -1 

♦ 0.015  (1  • SIN  (9- (1/2)  ) (4.7) 

He  then  have  as  the  target  cross  section  in  our  model  the 
following  expression; 

d = d x a (4.8) 

where  d and  U are  as  previously  shown. 

Fig.  6. a and  Fig  6.b  shows  an  'ideal'  pattern  and  a model 
pattern.  The  figures  given  by  Uriclc  are  for  1 yard  as 
reference  distance,  but  have  been  converted  to  1 meter 
reference  distance  in  Fig.  6.  To  go  from  dB(yard)  to 
dB (meter),  we  subtract  0.78  dB.  The  dB  as  given  in  this 
analysis  are  all  with  1 meter  as  reference  distance. 
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The  active  sonar  equation  is; 


2 

P -G  d G -X 
0 t - r 

P = Watts  (4.9) 

— z"1 — * 

(4TT)  -a 


P ■ power  received 

P * power  transmitted 
0 

G = gain  transmitting,  ref  Eg.  (4.1) 
t 

d = target  cross  section,  ref  Eg.  (4.9) 

G = gain  receiving,  ref  Eg.  (4.  1) 
r 

A = wavelength  in  meters. 

R = range  to  target  in  meters. 


This  may  be  rewritten  into  an  expression  of  power  received 
as  a function  of  the  variables  of  the  different  terms; 


SIN  (X  -If)  SIN  (X  -IT) 

t 2 2 2 2 r 

I I * a b c -R  | 

| | | J IT 

t r 


p = K 


2 2 2 2 T 

R • (a  -COS  t>  ♦ b -SIN  $) 


(4.10) 


K = the  product  of  all  the  constants  in  the  terms. 


For  more  detailed  discussion  about  gain, 
transmission  loss  and  reflection  (target  cross  section),  see 
[ 1 ; 110-1 11]  and  [ 6 ;29, 94, 263  ]. 

We  can  now  calculate  the  minimum  power  level  for  detection 
by  setting: 


33 


R^  = technical  detection  range 

X^  = 0 degree 

X =0  degree 
r 

<p  =90  degrees 
and  we  get 


P = K 
a in 


2 2 2 2 
a -b  -c  -R  *0 


4 2 2 

B -(b  ) 


(4.11) 


and  by  substituting  for  0 


P = K 
min 


2 2 2 2 
a -b  -c  -R  -3.08657 


4 2 2 

B *(b  ) 


(h. 12) 


a,  b and  c are  the  dimension  of  the  target  used  in  the 
model. 

We  assume  a 'standard*  target,  length  100  meters,  beam  15 
meters  and  draught  4 meters,  i.e. 
a = 100 
b = 15 
c = 4 . 


There  will  be  a detection  if  P/p  > 1.  Note 

min 

that  P/P  does  not  depend  on  K,  into  which  radiated  power 
min 

and  transducer  gains  have  been  included.  The  technical 


detection  range  R^  is  an  operationally  meaningful  surrogate 
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the 


for  these  parameters, 
"intensity  fraction". 


The  ratio  P/P  will  be 
min 


called 
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3 • Detection  Rule 


Any  intensity-fraction  calculated  during  a 
transmission  which  is  greater  than  1 is  a detection. 
However,  to  improve  the  model  at  close  ranges,  the  following 
modification  has  been  made  for  gain  variation  due  to 
relative  bearing. 

At  close  ranges,  the  relative  bearing  to  target  can 
alter  considerably  from  bow  to  stern.  Therefore,  the 
intensity  in  the  pulse  will  differ  along  the  target.  To 
average  this  intensity  both  for  the  radiated  pulse  and  for 
the  echo,  the  model  calculates  relative  bearing  to  the 
target  bow,  center  and  stern,  calculates  the  corresponding 
gain  factor  for  each  bearing,  and  finds  the  arithmetic  mean 
of  these  gain  factors.  These  two  average 
gainfactors  (transmitting  and  receiving)  are  then  used  in  the 
calculation  of  echo  intensity. 

The  model  does  not  recognize  a detection  unless  the 
tactical  situation  makes  it  possible  to  maintain  contact 
with  the  target  for  some  time.  To  be  precise,  the  following 
conditions  must  be  present; 

- torpedo  turn  rate  higher  than  bearing  rate 

- closing  speed  must  be  positive. 

- target  must  have  2 knots  doppler. 
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V.  PRESENTATION  OF  DATA 


A.  STOCHASTIC  ELEMENTS 

In  the  previous  description  of  the  model,  the  following 
input  values  are  stochastic; 
error  in  target  speed 
error  in  target  course 
error  in  target  range 

initial  torpedo  course  (not  aam  course). 


The  pnaary  stochastic  ir.  'he  torpedo 

performance  are  identified  as  error*  u target  speed  and 
course,  since  'h«  *•  . ' 'r.as'ic 

ones  used  in  coip  . t • - • » • : irst 

proble®  to  s«  »oi»«  « ..  erits 

in  order  to  reduce  » a . . • i • it  • . ie  as 

keeping  tae  resul* 

It  was  found  *aa'  a at  Iceized 

design  (raados  varia*-  • »...••»..►  ♦ :or. 

theprobabilityian.-  • • • .a.<ir*randoi 

variables,  using  * n-  » - . ♦ 

get  variates,  and  * hen  ri  to 

cover  all  cosbinations  : • • 


Sose  prelisinary  linlatlN  ru;;s  wer  . r-  :r  »hree 


versions;  complete  randomized  and  independent;  with 
antithetic  reduction  technique  (sectioning) ; and  the 
previously  described  procedure.  The  number  of  runs  needed 
in  order  to  keep  the  variance  low  for  the  result  was 
considerably  higher  for  the  first  two  versions. 
Accordingly,  we  selected  the  previously  described  procedure. 
It  was  found  that  a series  of  150  runs  was  sufficient  in 


order  to  give  a reasonable  accuracy  in  detection  probability 
and  at  the  same  time  keeping  the  total  CPU  time  for  a series 
of  runs  acceptably  low.  The  150  run  series  was  established 
by  dividing  the  range  of  probability  of  target  speed  errors 
into  15  equally  spaced  sections;  and  the  range  of 
probability  of  target  course  errors  into  10  equally  spaced 
sections.  Each  section  boundary  point  was  by  inverse 
probability  transformation  converted  into  a variate. 
Bearing  in  mind  that  speed  errors  are  normally  distributed 
and  course  errors  are  uniformly  distributed,  all 
combinations  of  target  data  error  are  plotted  in  Fig.  7. 
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J5.00 


Figure  7 - DISTRIBUTION  OF  ERROR  IN  TARGET  DATA 


Course  error  in  degrees 


B.  TYPE  OF  PRINTOUT  OF  DATA  AND  RESULT 


Each  series  of  150  runs  produces  a printout  as  shown. 

The  heading  cf  the  printout  gives  the  tactical  situation  and 
the  torpedo  parameters  in  the  given  run  series.  Also  , the 
printout  gives  the  sweep  lane,  which  is  the  width  of  the 
lane  where  the  torpedo  has  swept  through  by  its  sonar  lobe. 
The  coverage  ratio  gives  an  indication  of  the  fraction  of 
the  lobe,  which  is  covered  twice;  i.e.  how  much  the  lobe  is 
being  offset  from  its  previous  position  by  change  in  the 
torpedo  course.  The  question  of  offsetting  the  sonar  lobe, 
about  which  information  is  given  in  the  printout,  is 
discussed  later  in  Ch.  VI. 

Ref.  Fig.  8. 

For  each  run,  the  following  are  output;  target  data, 
torpedo  deflection  angle,  torpedo  main  course,  target  and 
torpedo  grid  position  at  end  of  run,  duration  of  torpedo  run 
and  length  of  torpedo  run. 

After  all  runs  in  a series  are  completed,  a summary  is 
given. 

Ref.  Fig.  9. 

The  summary  gives  detection  probability  for  a single 
detection,  2 successive  detections,  up  to  5 successive 
detections.  Also  mean  detection  range,  standard  deviation 
of  detection  range,  mean  aspect,  mean  detection  bearing 
relative  to  center  bearing  of  sonar  lobe  and  relative  to 
main  course  are  given. 


Lastly,  the  detection  range,  the  relative  bearing  to  the 
center  of  the  target  and  to  the  closest  part  of  the  target 
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at  detection  (relative  to  present  torpedo  course) , and  the 
target  aspect  at  detection  are  printed  for  each  run  for  a 
single  detection,  2 successive  and  3 successive  detections. 

It  also  should  be  noted  that  it  is  possible  to  get  a 
aore  detailed  printout  for  each  run  by  setting  IPEINT  = 0 in 
the  lain  program  (main  program  statement  035). 

Ref.  Appendix  D.  for  example  of  detailed  run  printout. 


Prom  the  printout  data,  it  is  possible  to  study 
different  aspects  of  the  detection  process  as  well  as  to 
generate  distributions  of  detection  range,  aspect,  bearing 
etc. 
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VI.  PARAMETRIC  TORPEDO  ANALYSIS 


A.  OBJECTIVES 

The  following  approach  was  used: 

The  torpedo  speed,  the  technical  detection  range  and  the 
lobe  width  were  assured  to  characterize  a torpedo  type. 
Within  the  type,  it  was  possible  to  change  the  turn  rate  and 
the  sweep  angle. 

A tactical  situation  was  characterized  by  the  attack 
angle,  the  target  speed  and  the  firing  range. 


The  following  questions  were  investigated: 

Can  a torpedo  be  improved  by  offsetting  its  sonar  lobe 
from  the  torpedo  heading  ? 

Rephrased;  it  may  be  asked,  is  the  sonar  lobe  searching 
in  the  right  direction  (most  likely  area)  by 
pointing  straight  ahead  along  the  torpedo  course  ? 

How  do  turn  rate  and  sweep  angle  affect  a torpedo's 
MOE  ? 

How  are  the  different  torpedo  types  related  to  each 
other  with  regard  to  detection  probability  (MOE)  ? 

In  the  analysis,  we  started  with  a reasonable  tactical 
situation;  target  speed  18  knots,  range  3000  meters, 
technical  detection  range  750  meters.  Initially,  we  changed 
the  attack  angles. 

With  regard  to  torpedoes,  we  started  with  three  types  of 
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torpedoes;  24  knots,  32  knots  and  40  knots;  all  with  20 
degree  lobe  width,  6 degree  per  second  turn  rate  and  30 
degree  sweep  angle. 

B.  OFFSETTING  SONAH  LOBE 

The  hypothesis  was  that  when  a torpedo  is  fired  on  a 
deflection  angle  course,  the  sonar  lobe  should  be  most 
effective  if  it  scans  across  the  bearing  to  the  target. 

Or,  the  sonar  lobe  should  be  offset  equal  to  deflection 
angle  (DA).  Ref.  Fig.  10. 


It  was  found  that  offsetting  had  a positive  effect  when 
attacking  from  ahead  of  target. 

Ref.  Fig.  11. a.  and  b. 

But  from  about  30  degree  to  about  110  degree  attack  angle 
the  effect  was  negative.  If  more  than  110  degree  attack 
angle,  there  was  no  effect. 

In  analyzing  the  fraction  of  offsetting,  we  analyzed  the 
case  of  30  degree  and  60  degree  attack  angle.  There  seemed 
to  be  no  effect  from  0.0  to  0.5  x DA;  if  more  than  0.5  x DA 
there  was  a decreasing  efficiency. 

This  was  found  for  2 types  of  torpedoes  (32  and  40  knots;  20 
degree  lobe  width)  at  2 different  sets  of  turn  rates  and 
sweep  angles. 


This  conclusion  applies  for  both  single  detection  and 
multiple  successive  detections;  however,  the  magnitude  of 
the  effect  is  changing  as  we  look  on  different  number  of 
successive  detections.  The  conclusion  was  that  there  is 
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little  to  be  gained  by  offsetting  the  sonar  lobe,  and  the 
sonar  lobe  was  therefore  not  offset  in  subsequent 
investigations. 


46 


Torpedo 
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EFFECT  OF  TURN  RATE 


r 


The  effect  of  turn  rate  was  investigated  in  the  range  3 
to  21  degrees  per  second  in  steps  of  3.  For  both  types  of 
torpedoes  the  model  showed  an  increase  in  MOE  as  turn  rate 
was  increased.  The  MOE  leveled  off  as  turn  rate  was 
approaching  15  - 20  degrees  per  second. 

The  reason  may  be  due  to  the  1 second  transmission 
interval  and  the  20  degree  lobe  width,  which  indicates  that 
the  torpedo  should  be  turned  at  a turn  rate  equal  to  lobe 
width  divided  by  transmission  interval  for  maximum  MOE. 
However,  as  the  number  of  successive  detections  required  is 
increased,  we  get  maximum  MOE  at  lower  turn  rates. 

Fig.  12  shows  the  change  in  MOE  with  turn  rate  for  30  and  60 
degrees  attack  angles. 

From  Fig.  15  where  different  combinations  of  turn  rates  and 
sweep  angles  are  plotted  versus  MOE,  we  see  that  the  effect 
is  negligible  from  about  60-80  degrees  to  180  degrees  attack 
angle. 


A 6 degrees  per  second  turn  rate  is  compared  with  what 
may  be  termed  an  'optimal'  turn  rate  in  Fig.  13.  The 
'optimal'  turn  rates  were  established  by  the  general  trend 
from  Fig.  12  and  Table  II. a and  II. b. 

The  following  turn  rates  were  identified  as  'optimal'; 

15  degrees  per  second  for  the  32  knots  torpedo 
18  degrees  per  second  for  the  40  knots  torpedo. 
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We  see  here  in  Fig.  13  a considerable  increase  in  MOE 
with  increase  in  turn  rate  for  attack  angles  less  than  60  - 
80  degrees  for  single  detection;  an  consistent  improvement 
for  2 successive  detections  in  the  same  area;  but  no  change 
or  a slight  detoriation  for  3 successive  detections. 

It  is  guite  obvious  that  a torpedo  which  requires  only  a 
single  detection  as  requirement  for  attack  has  a 
considerably  better  MOE,  and  a considerably  higher  potential 
for  improvement  by  changes  in  turn  rate,  than  a torpedo 
which  requires  more  successive  detections  for  classifying  a 
contact  as  a target. 
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D.  EFFECT  OF  SHEEP  ANGLE 

From  preliminary  simulation  runs,  it  was  found  that  from 
90  degrees  (inclusive)  to  180  degrees  attack  angle  the 
effect  of  the  sweep  angle  was  neglegible.  The  analysis  was 
therefore  done  from  20  to  50  degrees  sweep  angle  only  for  30 
and  60  degrees  attack  angle  for  both  the  32  and  the  40  knots 
torpedo. 

The  result  is  shown  in  Fig.  14. 


For  the  32  knot  torpedo  we  get  an  increase  from  30  to  40 
degrees  for  both  attack  angles.  From  40  to  50  degrees,  MOE 
either  levels  off  or  decrease  slowly.  As  a conclusion,  we 
established  40  degrees  sweep  angle  as  the  'optimal'  value. 
For  the  40  knot  torpedo,  the  MOE  was  fairly  steady  over  the 
whole  range  for  30  degrees  attack  angle.  For  60  degrees 
attack  angle,  there  was  a peak  at  30  degrees  sweep  angle, 
which  indicated  that  30  degrees  was  the  optimal  value. 


The  reason  for  the  different  sweep  angles  for  the  two 
torpedo  types  (Note;  both  have  6 degrees  per  second  turn 
rate)  may  be  due  to  the  time  it  takes  to  reach  the  target. 
The  shorter  time, the  less  area  on  each  side  of  the  main 
course  is  needed  to  be  covered  in  order  to  detect  a target; 
i.e.  a 40  knot  torpedo  needs  only  a 30  dergree  sweep  angle, 
a 32  knot  torpedo  needs  40  degree  sweep  angle. 
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The  reason  why  we  get  a peak  and  then  a reduction  in  HOE 
as  we  increase  sweep  angle  is  supposedly  due  to  a sharp 
decrease  in  speed  along  the  main  course  as  sweep  angle  is 
approaching  60  degrees. 

As  example,  for  a 40  knots  torpedo  the  model  gave  35  knots 
along  main  course  for  50  degrees  sweep  angle  as  compared 
with  38.6  knots  for  20  degrees  sweep  angle.  For  a slower 
torpedo,  the  effect  on  MOE  may  be  considerable  due  to  less 
speed  advantage  relative  to  the  target. 


E. 


EPFECT  OF  BOTH  SWEEP  ANGLE  AND  TORN  RATE 


In  the  previous  discussion  we  changed  either  sweep  angle 
or  turn  rate  for  both  types  of  torpedo  while  we  kept  the 
other  variables  constant. 


In  plotting  HOE  for  initial  torpedo  (32  knots)  value, 
optimal  value  for  sweep  angle,  optimal  value  for  turn  rate 
and  the  'optimal1  torpedo  (having  both  the  'optimal'  turn 
rate  and  sweep  angle),  we  get  Fig.  15. 

Observe  how  the  MOE  changes  as  we  apply  the  individual 
•optimal'  values,  and  the  MOE  obtained  by  applying  both  the 
'optimal'  values. 

At  this  stage,  no  trials  were  made  in  order  to  further 
increase  MOE  by  changing  sweep  angle  or  turn  rate  from  these 
values. 

One  essential  feature  is  that  virtually  none  of  the 
variables  so  far  have  had  any  effect  on  MOE  for  larger 
attack  angles  than  60-90  degrees. 


The  relative  difference  in  MOE  between  the  32  and  the  40 
knots  torpedo  types  is  shown  in  Fig.  16.  Both  torpedoes  are 
optimal  in  the  sense  that  the  best  values  for  turn  rate  and 
sweep  angle  have  been  choosen  for  that  specific  speed. 


Tactical  Situation 


Range  3000  m 
TA  Speed  1 S Knots 
Det.  range  750  m 


0°  30°  60°  90° 

120°  150°  18 

o 

One  detection  only 

Attack  angle 

□ 

Two  successive  detections 

A 

l hree » 

f TO  Speed  40  Knots 

32  Knots  1 

— < Turn  rate  1 8 °/s 

1 5 °/s 

l Sweep  angle  30° 

40°  J 

Figure  16  - COflPARISON  OF  TWO 

DIFFERENT  TORPEDOES 

omMMMP 


It  is  obvious  that  the  main  differences  are  for  large 
attack  angles;  more  than  60-80  degrees. 

Especially  if  the  acquisition  requirement  is  one  detection 
only,  however,  a 32  knots  torpedo  is  slightly  better  up  to 
60  degrees  attack  angle.  This  improved  MOE  for  the  slower 
torpedo  may  be  explained  by  a better  balance  between  the 
time  to  the  target  and  the  total  relative  speed.  A toe  high 
relative  speed  may  prohibit  the  torpedo  from  getting  the 
target  within  its  sonar  lobe  before  the  target  is  passed. 


Generally,  however,  the  higher  speed  torpedo  is 
superior,  especially  for  larger  attack  angles  (120  degrees 
and  more) ; this  can  be  explained  by  the  shorter  time  tc  the 
target. 


F.  EFFECT  OF  LOBE  WIDTH 


The  effect  of  changing  lobe  width  while  maintaining 
detection  range  is  shown  in  Fig.  17.  It  should  be  noted 
that  we  initially  started  the  simulation  with  an  'optimal' 
torpedo  with  20  degrees  lobe  width.  When  we  ran  the 
simulation  series  for  10  degrees  and  30  degrees  lobe  width, 
we  did  not  change  the  other  torpedo  parameters  in  order  to 
make  the  torpedo  'optimal'  for  the  new  lobe  width.  If  we 
had  carried  through  this  optimization  process,  we  might  have 
expected  an  increase  in  the  result  for  10  and  30  degrees 
lobe  width.  The  torpedo  parameter  in  question  would  most 
likely  be  tarn  rate,  ref  discussion  previously  on  page  61. 


The  interesting  points  from  Fig.  17  are; 

a 10  degrees  lobe  width  torpedo 
one-detection-only  acquisition  requirement 


with  a 
is  as  good 
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as  a 20  degrees  lobe  width  torpedo  with  a 
two-successice-detection  requirement.  This  should 
indicate  what  we  have  to  pay  in  additional  power 
transmitted  when  acquisition  requirement  is  high.  Or, 
where  to  invest  research  resources;  in  transducer  or  in 
echo  filtering. 

the  equally  shaped  curves  for  increasing  lobe  width. 
However,  we  also  observe  an  increasing  difference  in 
MOE  between  the  curves  as  attack  angle  is  decreasing, 
the  importance  of  the  correct  balance  between  turn 
rate  and  lobe  width  for  successive  detections.  He 
observe  for  a small  aspect  target  how  MOE  decreases 
drastically  when  we  reduce  lobe  width  from  20  degrees 
to  10  degrees  and  maintain  turn  rate  and  require  two 
successive  detections. 
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The  detection  range  is  a function  of  the  design  cf  the 
active  sonar  in  the  torpedo  as  well  as  sonar  condition  at 
the  time  of  the  torpedo  firing.  The  detection  range  as  a 
function  of  the  design  of  the  active  sonar  is  termed 
technical  detection  range.  The  detection  range  as  a 
function  of  both  the  design  and  the  sonar  conditions  is 
termed  tactical  detection  range,  or  just  detection  range. 
In  analyzing  the  detection  probability  as  a function  of 
detection  range,  we  assumed  optimal  sonar  conditions  by 
equal  technical  detection  range  with  detection  range. 

Detection  range  was  varied  in  discrete  steps:  375  - 750 
- 1125  - 150C  meters. 

Figs.  18. a.  and  b.  indicate  tnat  detection  probability 
is  a linear  function  of  the  detection  range  up  to  a 
detection  probability  of  0.8  -0.9  for  one  detection.  From 
the  model,  it  may  be  justifiable  to  approximate  the 
detection  probability  as  a linear  function  from  375  m to 
1125  m detection  range. 

From  the  model  and  the  given  assumptions,  thfr®  is 
little  usefulness  in  a homing  torpedo  with  less  than  300  m 
detection  range. 

The  same  situation  is  shown  in  Figs.  19. a.  and  t.  in 
another  cut  of  the  response  surface.  We  see  here  how 
consistantly  the  HOE  has  decreased  over  the  whole  range  of 
attack  angles  when  going  from  1500  m to  375  m detection 
range . 
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Otherwise  the  picture  in  Figs.  19. a.  and  b.  is  as  in 
previous  similar  figures;  a marked  decrease  in  MOE  with 
attack  angles  more  than  60  - 90  degrees,  and  with  the  faster 
torpedo  superior  over  most  of  the  range.  It  should, 
however,  be  noted  that  for  longer  detection  ranges  we  get 
maximum  MOE  at  0 degree  attack  angle  for  both  torpedo  types. 
This  effect  is  reduced  when  we  require  two  successive 
detections  for  acquisition. 

He  also  experienced  a considerable  decrease  in  MOE  for 
longer  detection  ranges  when  requiring  two  successive 
detections  instead  of  one.  It  seems  obvious  that  this 
reduction  is  due  to  a larger  lateral  movement  at  the  extreme 
range.  As  noted  previously,  we  increase  the  transmission 
interval  (increase  interval  in  order  to  allow  time  for  echo 
to  return)  when  the  detection  range  is  increased.  Keeping 
the  same  turn  rate,  the  sonar  lobe  will  turn  a larger  angle 
between  each  transmission,  which  can  have  a deteriorating 
effect  on  MOE  for  more  than  a one-detection-cnly  acquisition 
requirement. 
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Table  V - VARIATION  IN  DETECTION  RANGE 
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H. 


COMBINED  EFFECT  OF  LOBE  WIDTH  AND  DETECTION  RANGE 


The-  following  approximate  relationships  exist  between 
lobe  width,  detection  range  and  sonar  power: 

v 2 

P G <3  G -A 
0 t r 

p = Watts  (4.9) 

3 — H 

(4  1 T)  -R 

where 

G = g = (4  •IT  /w)  (6.  1) 

t r 

2 

w = L 

L = 2 x lobe  width. 

Ref  [ 1 ; 4 9 ]. 

w is  defined  as  solid  angle.  The  given  equation  is  valid 
for  small  lobe  width  only.  For  larger  lobe  width  the  exact 
relationship  is; 

w = 2 *1T  • ( 1 - COS  1)  (6.2) 

1 = lobe  width. 

The  approximate  relationship  is  close  enough  up  to  60 
degrees  lobe  width. 

By  substituting  the  approximate  relationship  into  Eg. 

4 

4.9  , we  get  a reduction  of  L in  receiving  echo  due  to 

change  in  lobe  width, 
or 
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4 


constant. 


(6.3) 


(L-R)  = 

which  coBDine  range  and  lobe  width,  and  implies  that 
detection  range  is  inverse  proportional  to  lobe  width  for 
constant  power  transmitted. 

It  is  therefore  possible  to  plot  this  function  for  constant 
power  transmitted,  and  use  this  as  a prediction  of  hew  MOE 
may  change  with  change  in  these  two  torpedo  parameters  (lobe 
width  and  detection  range)  . 

This  is  done  in  Fig.  20;  and  indicated  by  the  dashed  line 
going  through  20  degrees  lobe  width  and  750  m detetion 
range. 

We  then  ran  some  simulation  series  in  order  to  generate 
data  points  from  the  model.  The  data  points  gave  the  MOE  , 
and  by  fitting  curves  we  were  able  to  get  some  indication  of 
the  relationship  between  the  lobe  width  and  the  detection 
range  as  given  by  the  model. 

The  application  could  be  as  follows; 

For  a given  torpedo  with  lobe  width  20  degrees  and  a 
detection  range  of  750  m,  we  ask  the  question,  can  MCE  be 
increased  without  increasing  power  transmitted  ? 

The  dashed  curve  through  the  point  (20  degrees,  750  m)  is  a 
constant  power  curve,  and  by  following  the  curve  we  observe 
how  MOE  is  changing. 

From  the  figure,  it  is  obvious  that  a narrower  lobe  and  a 
longer  detection  range  gives  a better  result.  But  we  also 
observe  the  assymptotical  feature  of  the  curves.  We  reach  a 
point  where  the  constant  power  curve  and  the  constant  MOE 
curve  are  parallel. 

However,  it  should  be  born  in  mind  that  the  theoretical 
relationship  between  lobe  width  and  detection  range  is  an 
approximation  which  does  not  account  for  absorption-effect 
or  surface-effect.  This  implies  that  the  constant  power 
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curve  in  real  life  will  be  lowered.  Only  a more  detailed 
analysis  can  say  how  much. 


I. 


EFFECT  OF  FIRING  RANGE 


The  most  important  factor  in  achieving  high  detection 
probability  is  the  difference  between  estimated  target 
position  and  actual  target  position  at  the  time  when  the 
torpedo  is  in  position  to  detect.  The  effect  on  the 
detection  probability  is  mainly  due  to  the  time  the  torpedo 
takes  to  reach  within  detection  range  of  target  and  the 
speed/course  errors  in  target  data. 

As  we  increased  the  firing  range,  we  experienced  as 
anticipated  a degradation  in  MOE.  This  degradation  was 
experienced  for  both  the  32  and  the  40  knots  torpedo. 

The  variation  in  firing  ranges  were  at  the  following  values: 
1500  - 3000  - 5000  - 7000  meters. 

Fig.  21. a.  and  b.  shows  consistently  the  importance  of 
short  firing  ranges.  This  applies  to  both  one  detection  and 
two  successive  detections. 

Fig.  22. a.  and  b.  shows  an  additional  advantage  with  short 
firing  ranges;  a considerable  improvement  at  firing  with 
small  aspect  (attack  angle),  less  than  30  degrees.  Again 
this  applies  for  both  types  of  torpedoes. 

Also,  we  get  an  indication  that  at  short  ranges,  about  1500 
meters,  there  is  no  significant  difference  in  MOE  of  the  two 
types  of  torpedoes  up  to  an  attack  angle  of  90  degrees. 


82 


Tacticai  Situation:  Torpedo  Parameters: 

TA  Speed  18  Knots  TO  Speed  32  Knots 


1500  3000  5000  7000  m 

Firing  range 


Figure  21  - EFFECT  OF  FIFING  RANGE 
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Detection  Probability 
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Tactical  Situation  Torpedo  Parameters: 
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Torpedo  Parameters  Tactical  Situation 


TO  Speed 

32 

40  Knots 

TA  Speed  18  Knots 

Lode  width 

20 

20  0 

Det.  range  750  m 

Sweep  angle 

40 

30  0 

Turn  rate 

15 

18  °/s 

Two  successive 
detections 


Figure  22. b.  - COMPARISION  OF  TWO  TORPEDOES  WITH  CHANGE 

IN  FIRING  RANGE 
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Tactical  situation  Torpedo  parameters 
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87 


J.  EFFECT  OF  TARGET  SPEED 

Generally,  we  anticipated  a degradation  in  MOE  as  the 
target  speed  was  increased.  And  overall,  this  was 
confirmed. 

The  simulations  were  carried  through  at  12,  18,  24,  30  knots 
target  speed. 

However,  fig.  23. a and  b shows  some  interesting  patterns 
regarding  optimal  attack  angle  for  different  target  speeds. 
For  a 32  knots  torpedo,  at  60  degrees  attack  angle,  the 
torpedo  is  equally  good  for  any  type  of  target  speed  for  one 
detection  only.  For  two  successive  detections,  the  torpedo 
is  equally  good  between  30  and  90  degrees  for  12  and  18 
knots  target.  A 24  knots  target  gives  a consistently  lower 
MOE  over  the  whole  range  of  attack  angles,  and  the  2 
simulation  runs  with  a 30  knots  target  confirmed  that  trend 
for  the  32  knots  torpedo. 

We  may  form  the  conclusion  that  for  one  detection  only 
60  degrees  attack  angle  is  an  optimal  attack  angle  for  the 
range  of  target  speeds.  For  two  successive  detections,  30 
to  90  degrees  attack  angle  gives  equally  good  MOE  between 
target  speed  of  12  and  18  knots. 

One  interesting  point  is  that  it  seems  that  if  the 
target  speed  is  less  than  0.4  of  the  torpedo  speed  the 
optimal  attack  angle  shifts  forward  to  0 degree. 

This  also  applies  to  two  successive  detections. 
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Tactical  Situation: 
Range  3000  m 
Det.  range  750  m 


Tactical  Situation: 
Range  3000  m 
Det.  range  750  m 


Torpedo  Parameters: 

TO  Speed  40  Knots 
Lobe  width  20° 
Sweep  angle  30° 

Turn  rate  187s 


One  detection 
Two  successive 
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For  a 40  knot  torpedo,  in  addition  to  the  point  of 
optimal  attack  angle  at  0 degree  for  slow  target  speeds,  we 
also  experienced  a relatively  low  MOE  for  slow  targets  in 
the  range  45  to  105  degrees  attack  angle,  compared  to  fast 
targets.  But  as  a compensation,  MOE  is  increased  for  small 
attack  angles  and  the  astern  attack  angle  compared  to  fast 
target.  Obviously,  some  type  of  a breaking  point  is 
experienced  for  target  speed  of  .4  or  less  of  torpedo  speed. 

Why  a slow  target  produces  this  increase  in  MOE  in  the 
two  extreme  cases  (ahead  and  astern)  may  be  explained  by  the 
balance  between  time  to  reach  detection  range  and  the  total 
relative  speed.  It  is,  however,  more  difficult  to  give  any 
explanation  of  why  a slow  target  should  produce  a lower  MOE 
for  some  attack  angles  than  a faster  target  does.  One  wculd 
have  anticipated  an  increase  in  MOE  over  the  whole  range  of 
attack  angles  for  a slow  target. 
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Detection  Probability 


Torpdo  Parameters:  Tactical  Situation: 


TO  Speed 

32 

40  Knots 

Range 

3000m 

Lobe  width 

20 

20° 

Det.  range 

750m 

Sweep  angle 

40 

30° 

Turn  rate 

15 

187s 

Figure  24  - COMPARISION  OF  TWO  TORPEDOES  WITH  CHANGE  IN 

TARGET  SPEED 
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VII.  TACTICAL  ANALYSIS 


In  addition  to  the  detailed  parametric  analysis,  which 
has  been  shown  previously,  we  also  could  expand  the  analysis 
to  cover  more  tactical  related  problems.  If  we  assume  a 
given  target  speed,  we  could  construct  detection  probability 
charts  as  shewn  in  Fig.  25.  a.  and  b. 

This  analysis  would  then  naturally  fall  into  two  areas: 
direct  comparison  of  two  or  more  different,  types 
of  torpedoes. 

effect  of  tactical  situation  on  the  detection 
probability. 

The  two  charts  (Figs.  25. a.  and  b.)  were  formed  by 
running  simulation  runs  for  differnt  tactical  situations 
(range  and  attack  angle) , and  then  fitting  constant 
detection  probability  curves  through  the  data  points. 

The  use  of  these  types  of  charts  falls  into  two  areas: 
Evaluate  different  torpedo  types  for  different  tactical 
situations;  essentially,  which  torpedo  is  best.  Or  for  a 
given  tactical  situation,  how  could  the  situation  be 
improved,  and  what  options  exist. 

The  first  type  of  use  applies  mainly  to  operational 
planning;  operational  requirement  in  the  design  phase  of  a 
torpedo  and  procurement.  By  laying  one  chart  atop  of  the 
other;  we  get  a visual  picture  of  how  much  is  improved  when 
using  a 'better*  torpedo,  and  for  which  tactical  situation. 
The  shaded  area  in  Fig.  25. b.  shows  how  many  more  tactical 
situations  have  been  covered  when  going  from  a 32  knots 
torpedo  to  a 40  knots  torpedo  for  0.25  detection 
probability. 
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Det  range  750  m TO  Speed  40  Kno 

TA  Speed  18  Knots  Sweep  angle  30° 

LoPe  width  20  ° 
Turn  rate  1 8 


Figure  25. b. 
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The  other  type  of  use  of  the  charts  is  tactical.  When  a 
firing  unit  decides  to  attack,  and  finds  itself  in  a given 
tactical  situation,  the  question  is:  What  to  do  ? 

For  given  target  speed  and  own  max  speed,  the  charts 
make  it  possible  in  a simple  way  to  decide  where  to  go  and 
what  course  to  keep.  Also  from  the  charts,  one  can  decide 
where  on  the  relative  course  is  the  optimal  firing  position. 
For  a submarine  attacking  a zig-zagging  target,  the 
Commanding  Officer  can  better  make  his  evaluation  of  when  to 
fire,  as  the  attack  angle  and  the  distance  are  continuously 
changing.  He  can  see  what  improvement  to  expect  when  the 
target  will  change  course  next  time.  An  example  of  a 
tactical  situation  and  the  course  of  action  to  follow  are 
given  in  Fig.  25. a. 

These  points  also  bring  up  the  question  of  what  to 
improve  in  the  operational  picture;  the  firing  unit's 
ability  to  acheive  a good  firing  position  or  the  torpedo's 
ability  to  detect  target  from  non-optimai  situations. 

In  this  discussion,  the  guided  torpedo  has  to  be  brought 
into  the  picture.  The  effectiveness  of  guidance  has  not 
been  adressed  at  all  in  this  study,  basically  because  that 
would  have  significantly  expanded  the  scope  of  the  study,  as 
well  as  bringing  in  the  whole  problem  of  fire  control 
equipment,  its  effectiveness  and  its  reliability. 
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VIII.  CONCLUSIONS 


The  study  was  carried  out  in  order  to  investigate  the 
detection  process  of  an  active  sonar  homing  torpedo  used 
against  surface  ships. 

Specifically,  we  wanted  to  study  the  effect  of  changes 
in  torpedo  parameters  such  as  torpedo  speed,  turn  rate, 
sweep  angle  and  detection  range,  as  well  as  changes  in  the 
tactical  situation  such  as  target  speed,  firing  range  and 
attack  angle. 

In  an  attempt  to  gain  insight  into  the  complexity  of  a 
homing  torpedo,  the  described  model  was  built  and  the 
simulations  done  as  previously  shown. 

In  designing  a homing  torpedo  and  evaluating  torpedo 
tactics  the  detection  probability  is  an  essential  part  of 
the  total  effectiveness  of  the  torpedo. 

To  be  able  to  hit  the  target,  the  torpedo  has  first  to 
detect  it,  which  justifies  why  we  started  out  with  analyzing 
the  detection  process. 

Also,  as  part  of  this  analysis  we  investigated  certain 
aspect  of  the  next  step  in  the  operational  process; 
acquisition. 

It  is  not  difficult  to  visualize  tests  which  may  be  used  in 
order  to  recognize  an  echo  as  a detection  and  subsequently  a 
target  to  attack.  Some  of  these  tests  may  be  doppler, 
successive  detections,  detections  within  a given  range, 
2-of-3  detections,  size  of  echo,  length  of  echo  etc. 

The  problem  of  false  echo,  however,  was  net  approached  in 
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this  study.  That  would  have  to  be  the  next  to  consider  in 
relation  to  reducing  the  number  of  successive  detections  in 
order  to  acquire  a target. 

In  order  to  allow  for  the  errors  in  tracking  of  the 
target  before  firing  and  also  small  maneuvering  of  the 
target  after  firing,  we  introduced  errors  in  the  target 
speed  and  course  when  calculation  of  the  torpedo  main  firing 
course  was  done. 

During  runs  the  torpedo  was  unguided,  and  did  not  react 
on  any  detection;  i.e.  it  did  not  attack  the  target.  For 
sonar  condition,  isovelocity  was  assumed  and  no  surface 
effect  was  built  into  the  model. 

The  result  gave  certain  insight  into  the  complexity  of 
the  detection  process,  stressing  the  importance  of  a good 
tactical  firing  position,  and  ' high  speed  torpedo  with 
long  detection  range. 

However,  the  data  also  showed  the  relationship  between 
detection  range,  lobe  width  and  turn  rate,  as  well  as 
weighting  the  sweep  angle  in  relation  to  torpedo  speed.  It 
can  also  be  concluded  from  the  results  that  changes  in 
torpedo  parameters  as  turn  rate  and  sweep  angle,  which  may 
be  inexpensive  modifications,  will  not  give  a significant 
improvement. 

Generally,  the  overall  important  factor  was  the  time  the 
torpedo  used  in  order  to  travel  within  the  detection  range 
of  the  target.  This  was  due  primarily  to  the  error 
generated  in  the  target  data,  obviously  the  actual  value  of 
the  result  is  sensitive  to  these  assumptions. 

However,  the  understanding  and  insight  in  the  detection 
process  achieved  by  simulation  should  not  be  reduced  by 
other  assumptions  with  regard  to  error  in  target  data. 
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with  regard  to  the  analysis,  the  result  has  shown  a 
consistent  and  general  trend  that  if  we  are  able  to  require 
only  one  detection  for  acquiring  a target,  the  detection 
probability  is  significantly  higher  than  if  more  than  one 
detection  is  required.  And  what  is  more  important,  the 
potential  for  improving/optimizing  a homing  torpedo  is  also 
significantly  higher  for  one  detection  only. 

This  implies  a large  payoff  for  other  methods  of  keeping 
down  the  probability  of  false  detections. 

Secondly,  a high  speed  torpedo  has  shown  a general 
superiority  in  MOE.  This  was  specially  obvious  in  attack 
angles  greater  than  90  degrees,  which  tends  to  make  a high 
speed  torpedo  more  of  an  all-round/reliable  torpedo  with 
regard  to  tactical  situations. 

Thirdly,  except  for  changes  in  attack  angle  and  firing 
range,  the  detection  range  seems  to  influence  the  HOE 
strongly. 

These  three  remarks  all  point  towards  an  improvement  in 
the  sonar-/filter ing-area  as  the  most  promising  area  in 
which  to  carry  out  research  and  invest  effort. 

This  study  has  also  pointed  out  the  advantage  of  high 
torpedo  speed  and  firing  at  short  ranges.  There  exists 
therefore  considerable  argument  for  a short  range,  high 
speed  torpedo,  given  that  one  is  able  to  position  the 
torpedo  at  a short  firing  range;  i.e.  a small,  simple 
torpedo. 

Basically,  there  are  two  schools  of  thought; 

a highly  sophisticated  torpedo;  long  range,  guidance, 
expensive,  but  close  to  the  one  shot-one  hit  idea, 
a simple,  high  speed  torpedo;  short  range. 
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non-guidance,  inexpensive,  and  requiring  either  a 
firing  unit  which  can  get  into  an  optimal  firing 
position  or  a larger  number  of  shots  to  acheive  hit. 


The  result  may  be  useful  in  giving  example  of  how 
tactical  guidelines  can  be  evaluated  by  the  simulation 
approach.  But  more  significant  is  pointing  out  the 
importance  of  torpedo  capability  and  the  tactical  situation. 
Obviously,  we  have  to  look  on  the  whole  torpedo  system, 
including  the  firing  unit.  Investment  in  resources  and 
effort  should  not  necessarily  be  spent  only  on  the  torpedo 
in  order  to  increase  its  effectiveness,  but  may  be  spent  on 
the  firing  unit  as  well  in  order  to  make  the  unit  able  to 
reach  a tetter  firing  position. 

A follow-on  of  this  study  may  be  to  investigate  the 
attack  process  of  the  torpedo,  including  the  acquisition- 
and  hit- problem. 

Then  the  question  of  guidance  during  torpedo  run  should  be 
analyzed  in  order  to  better  evaluate  the  problem  of  choosing 
between  a few  sophisticated,  expensive,  guided  torpedo 
system  or  many  simple,  inexpensive,  nonguided  torpedo 
systems. 
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APPENDIX  B 


FLOW  CHART  FOR  SIMULATION  PROGRAM 


A TOE  EEOC  SIMULATION. 


MAIN  EBOGEAM. 


PACE  1 


A IC&PECO  SIMULATION. 


A 10 B P £ CO  S ISOLATION . 


MAI N PBOGBAM. 


SIMULATING  AN  AC1IVE  dOfll  NG  TORPEDO  CUBING  SEASON. 

THE  PBCG8AM  IS  BUN  IN  0.5  SEC  STEPS. 

COMMON  ISEEDi.  TTIME,  TO.  TA  , T3ATE,  BANGS,  ALEA,  LAMB  C.  TAOEC, 
oiAfi,  BAC,  TAC.  CCCB,  CEVSP,  ENG,  EN,  PH2,  MCCUHS,  TCOOBS, 

MET,  MXM.  IK.  IDTIMS,  ITIM2,  XT,  YT,  XIAB,  TTA3,  TJIST,  MOIST, 
TUENTO,  INTVAL,  PHI,  B MAX  (5,6) , CHANGE,  OlSl.IPBINT 

COM MON/O A I A/ I A , OPLGB 

CCMMON^TABGET/TACMG,IAM1 .BNGMOD ,OA1,COUR,CE (10)  ,SE(15)  ,JBUN, 

DIMENSION  SM  (5,5)  , VAB<5),  OET(150,S),  DETB  ( 150,5)  , STO  (5)  , 

ASEEC  (150,5) ,d£BB (150, 5) .CLOSE (1 50,5)  , KON ( ISO, 5) 

BEAL  LAMEC,  MCOUBS,  MXT,  MXH , MOIST,  LAMBDG 

I SIEGES  SONOUT 


SETTING  Of  CONSTANTS  (STEP  1) 


CALL  OVPLCU. 

__J 

E HI  -3.14159265* 
EH*  -2. -PHI 
BAD  -EH2/360. 
ISEE01  -362716 
I5EED2  -S616S5 


CCOBEC  - MAX  SBBOB  IN  TABGET  COUESE  ESTIMATE 


CCOaEC  -15. 

CCOB  -CCOE  EC-BAD 


CSPEEO  - ST  0 E 1 IN  TARGET  SPEED  ESTIMATE 


CSPEED  *3. 
CEVSE-CSEEEC-0.5 
ICTIME  -1 


SET  NOMEEB  Of  IliBATIONS 


| IBUN  -150 


(CONTINUED  ON  E 


Li  21 
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A TOEPEDO  SIHOlAIIOtl.  3 AIB  P80GBA3.  PAGE  2 

I 

SET  PBI3T  OUT  3CDE 

l 

1 I PH  1ST  »T I 

I 

SET  LOSE  0?F  TOE  PE  CO  CENIEB  3EABING 

1 

I CfLOo^O.  ] 

( 


SET  T1EI.ES  TC  ZEEO 

1 

(STEP  2) 

s' 

♦ 

♦♦♦♦♦♦♦♦♦♦  11 

♦ 

♦ ♦ 1*1, 5 

♦ 

♦ | 

♦ 

♦ ♦ DC 

♦ 

♦ 12 

♦ 

♦ + J»1,5 

♦ 

♦ 1 

♦ 1 

♦ 1 EH(I,J)  =0.  1 

1 1 

♦ 

♦♦♦♦♦♦♦ +CCS1INUE 
♦ 

♦♦♦+++++COBTIBOE 

1 

♦ DC 

♦ 

♦+♦♦♦♦♦♦♦♦  13 

♦ 

♦ ♦ J«1,S 

♦ 

♦ 1 

♦ 1 

♦ * CO 

♦ 

♦ ♦♦♦♦♦♦♦♦♦  14 

♦ 

♦ ♦ I=1,IJUN 

•f 

+ 1 

♦ 

♦ I XC3(I,J)  *0 

1 

♦ 1 

♦ 

♦♦♦•*  ♦♦♦♦continue 

♦ 

♦♦♦♦♦♦♦♦CONTINUE 

1 

| LA2  >0.  1 

(CONTINUED  ON  r 1g £ 

3) 
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1 TOHPEDO  SIMULATION. 


MAIN  p HOG  BAM  . 


PAGE 


COMPOTE  IlHGET  E a SOBS  AND  STOBE 

i 


♦ DO 

► 20 

♦ I»1. 10 


| CE  (I)  *-  (15.  *HAD*CC03/1 0 . ) ♦ (CCOS 

/ 10. ) *1*2. 


20  ♦♦♦♦♦♦♦♦CCSIINUE 


SE  ( 1 

*- 1. 8 7 37*DE  7SP 

i SE 

2 

*- 1. 2825*DE?S? 

1 SE 

3 

=-C . 9 6 6 *DE  V S? 

SE 

4 

»-C. 72S5*DEVS? 

SE 

5 

*-C.52S*DEV5P 

S E 

6 

*-0. 34«CEVSP 

SE 

7 

»-C.  1e75*DEVSP 

SE 

8 

•0*DE1iS? 

SE 

9 

= 0. 1 6 7 £ *DE V SP 

SE 

1 

=0. 34»DE7S? 

SE 

11 

»0. £2£»DEVSP 

SE 

12 

»0. 7285«DE7S? 

SE 

13 

=0 . S6fi*DE VSP 

1 SE 

14 

= 1. 2S2S*D2VSP 

1 SE 

15 

*1.87  37*DEV5? 

100 


♦ DO  ♦ 

♦ 900  ♦ 

♦ ♦ JBUN*1,I£ON  ♦ 

I 

♦ SET  SON  COUNTESS  (STEP  3) 


IFLAG=1;  TOC  LOM  ICEP230  SPEED 

I 


IFLAG-0 
ICC  NT- 0 
JCONT-O 
III 3E=  0 
IL  *C 
IK  *0 
CIST  =0. 
311  »0. 

3X3  *0. 

SUNOUT  =0 


(CONTINUED  ON  P iG  1 


la  i 


4) 
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I TOBPEDC  SiaOLAIICN. 

♦ 

♦ 


ana  pbogbab. 


PAGE  4 


HEAL  IN  SETTING (TOBP  AND  TACTICAL) 


FISSI  SUN  ? (STEP  4) 


IF  * 

J80N.GI.1 


I 160  | 


CALL  PAH  NET 


•••SHITE (6,128) OF LOB 


?CS»»I(1X,/,1X,'SCNA-B  2AIN  LC3Z  OPF-SET  FBC8  CZNTS3  ' 
' ESAHING ',16.2,'  II aSS  DEFLECTION  ANGLE',/) 

| 

i IVISI=IC/2. 

I 2DIST=IA/2. 

I IUBNTC  =TB  ATE/2 . 

INTVAL  = IFIX  ( (TTII1E/C.5)  *0.5) 


INTVAI  GIVES  NO a EES  OP  TiaSSTSPS  FOa  EACH  TBANS8ISSI0N 


. » * . 

■ IF  * 

IPBISI.EC.O 


PBINT  CF  HEADING 


•••MBITS (6, 220) 


F0EaAT(1X,//,1X,<aUN' ,4X-'EST  OP  T ARG ST* , 7 X , 1 TOHP  TCBP  8', 
SX.'TCHP  CCC3C,5X, 'TAHGEr  CCOED  SUN' , 32 , • TOBP • , / , 

1 X f 1 NO  COOSSE  SPcED  BANGS  DA • , 4 X, 'COOaS E* , 1 0X , ' X • , 6 X, 
'X1 ,7X,' X',6X, '!• ,6X,' STOP' ,2X, 'BUN' ) 

(CCNTINUED  ON  picE  5) 
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A TOHPECO  SIMULATION  • 
♦ 

♦ 

♦ 


MAIN  PHOG8AM . 


160 


CALCULATE  TOBPECC  START  POSITION  (STEP  5) 

I 


2TAB  *1S00C. 
HAS  =1500C. 

C P A 1 = EANG E 
£ MG  »IAC*ESAB 


I 

. * * . 

* IE 
ENG.GI.PB2 


. * 


| ENG=BNG-PH2 


. • * . 

• IE  » . 

3NG.L1.C.  • 

. * 


• . * 
!— 


3N=PH2*BNG 


II 

XT 


«XTAS  + aAMGE»Si:i  (3MG) 
•ITAB+ RANGE “COS  (3NG 


198  ♦ 
♦ 
♦ 
♦ 
♦ 


152  ♦ 


t xp  < 

IP8IS1. EC. 1 


• ••MBITS  (6,  158)  JBON 

?0HMAT(1I,/,6l,'BaN  NUMBER 


I 152  | 


.14) 


CALCULATE  TOE  EE  TO  DEFLECTION  ANGLE  AND  FIB  IN  G SITUATION 


CALL  RISING 
(CONTINUED  ON  pIgE 


6) 


PAGE  5 


12  U 


♦ ♦♦♦♦♦♦ 


A TOEPEDO  SIHOLAIION. 


sain  progbas. 


SET  DETECTICNTASIES  TO  ZEaC  (STEP  6) 

♦ 

♦ I 

♦ 4 ~4 

♦ ♦ DC  ♦ 

♦ 10  4 

4 4 1*1,5  ♦ 

♦ I 

♦ 

4 4 DO  *• 

44444444+4  9 ♦ 

4 4 J3  1*6  * 


l caAZ(i,J)  -o. 
) 


9 +4++4+4+CCSIXH0E 

4 

10  44++44+4CC5TIHUE 


TEST  TOEPEDOaCS  CCT  (STEP  7) 

I 


( THUD  *TKANGE/TO 
J CIOS?-IA»CCS (SEAS) 4T0»C0S (DA) 
| BBUM  »flASGE/CLOS? 


4 . * • . 

4 . • If  * . 

4 * ISDN. GT. HBOS  • 

4 * . . 

+ * . .4  T I 499  | 

4 * . * 

1 

4 I 

4 ***«aiI2(6,;i1) 

231  4 PCBHAT(1X,/, IX. ’TARGET  OOTSIDE  TOBPEDO  BANGE’) 

4 

4 

4 

♦ Z2SC  TABLES 

4 

♦ I 

4 * 4 

4 4 DO  4 

444+44444+  23  ♦ 

4 4 0*1,5  ♦ 

4 4-— — — 4 

4 I 

4 | CITE  (JBCN  , 0).  *0^  | 

(CONTINUED  ON  tld  7) 


PAGE  6 
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1 TORPEDO  SlfiULATICH. 
♦ 


BAIN  E HOG BAH . 


PAGE  3 


i 

C All  DETECT 


CHECK  CPA  (CLOSEST  POI3I  OP  APPROACH)  (STEP  8) 


I 


I IL 


«U*1 


I 


II 

IL. 12. 20 


* „ 


I 500  | 


J 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 


♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ ♦ 
♦ 


IL  * C 
CPA  *CIST 


• * * • 

. * II  * 

* CPA 1-CFA. L2.0. 

* . 

• . . * 


l a(JNOU?*i 


| CP A 1 *Cf A 


I SCO  | 


990 


GEHEHA12  STA IISTIC S (STEP  9) 
CCBIXMUE 

I 


♦ * DO 

♦♦♦♦♦♦♦♦  910 

♦ IKl»t,5 


I 

(CCSTIM02D  OH  ploE  9) 
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A TORPEDO  SIEOLAIICN.  NAIN  FaOGBAd. 

♦ 

♦ I 

♦ 

♦ 

♦ 

♦ KON<  ) « DETECTIC3/NO  DETECTION 

♦ 

: [ 

♦ . * • . 

♦ . • ii  * . 

♦ * sen  (ixi.i)  . ct.  i.  * 

t » . . * r | KOB  (JBUS,IKl)»1 

♦ • . • 

: * I 

♦ 

♦ 

♦ 

♦ DET  ( ) * DISTANCE  TO  TARGET  AT  DETECTION 

♦ 

♦ I 

♦ | LET  (JHGM  ,IKI)  » B 2JUC  ( ZKL , 1 ) I 

♦ — — — 

♦ I 

♦ 

♦ 

♦ 

♦ OETE  ( J * SEARING  TO  TARGET  AT  DETECTION 

♦ 

♦ \ 

♦ 1 DET3(JBDN,IKL)  »BMAXjIKLr2)  | 

♦ I 

♦ 

♦ 

♦ 

♦ 

♦ ASPEC  ( ) » TARGET  ASPECT  AT  DETECTION 

♦ 

♦ ASPEC  / JR  CB ,1  XL)  «BNAX4IKL#5) 

♦ 

♦ 

♦ » 

♦ CLCS3<  ) * BEARING  TO  CLOSEST  PABT  0?  TARGET 

♦ 

♦ CLCS3  (jaOtf,IAL) *B  a AX (IKL,  4) 

♦ 

♦ 

♦ 

♦ DEBE  ( ) * BEL  EE  A BING  PROS  .IAIN  TOBP  COUBSE  TO  1ABGET 

♦ 

♦ J 

♦ | DEBE (JBON,IKL)  *B«AXJIXL^6) | 

♦ I 
910  ♦♦♦♦♦♦♦♦CCATINOE 

(CCBTISDED  OB  lid  10) 


? AG3 
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PASS  10 


A ICS  PECO  SIBOLAXICN. 
♦ 


BAIN  PBOGEAB. 


. • 

• IP 
IPLAG.EQ. 1 


I 991  | 


PBIMI  COT  (SUP  TO) 

I 


| COOS  *BCOOFS/RAD 
TC  -ICCCES/RAD 
IT  »IIiaE/2 


• IP  * 

IPBI XT  . EC • 0 


I 195  ( 


232  : 


•••WHITE  (6  , 232)  JR  ON  ,T  AC*1G  / rAB1#3NGBCO/OA1,COOH,XT#TT#  IIAa.ITAH, 
IT  , BIT 

PCSaATpX./.II.n.OX.PS. 1t2X,F4. 1, 1X,F6.0,2X,F5.1,3X, 

P5.  1,?X,Pt.fi,  1X,Pc.O.,2X,?o.0, 1X,Fo.0,3X,li,2i,F6.0) 


. • IP  • 
• EHAX  (1  , 1)  . PC.O. 


991  | 


195  ♦ 

200  ♦ 
♦ 


| 9C0  | 

•••WHITE  (6,200)  IT 

FCEHAI (1 1. *aOS  STOPPED  APTEE  '.IS,1  SECONDS',//, 
IX, 'BUN  DATA  AS  FOLLOWS  IT  END  0?  30.N') 

(CCNTI.10ED  ON  pIgS  11) 
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A TORPEDO  SiaULAIION. 
♦ 


3AIN  PROG  BA. 1 . 


PAGE  ' 11 


♦ 

♦ 

♦ 

202  t 
♦ 
♦ 

204  ♦ 
♦ 
♦ 
♦ 
♦ 


•••WHITE (6, 202) HX1.0IST 

FORBATflX,' TOTAL  TORP  BUS  • , P9 . 1,/ , 1 1 , • DIST  TO  IABGES  ',F9.1) 

***88112  (6, 2C4)  XT.YT,XTAa,YTAR,COUa,TC 

FOHBAT(1X,'TOSP  X-COORD  ' , 2X ,?9 . 1 , 4X , ' TOHP  Y-COORO  *,2X,F9.1,/, 
IX, 'TARGET  X-CCORD  ', F9. 1 , 4 X TARGET  Y-COORD  '.F9.1,/, 

11,‘TOaP  SAIH  COOBSE  *,F9.3,eX, 'TOBP  COURSE  »,1a,F9.3) 


. * IE  * 
HHAX  (1,  1)  .EQ.O. 


I 991  l 


*»*BEITE  (6,206) 

•••WRITE  (6,203)  (I,  (RKAX,(I,J)  ,J*1,5)  ,1*1,5) 


’CiSTES1 ,5X, '3ASS2  - CLJS2ST1 ,4X, • ClOSSST' , 
8X,’ ASPECT') 

FORMAT  <3X,I4,7X,Pie. 1 , 7X , ?8 . 2 , 7X, F 10. 1,9X. 
P8.2,7x,E?.2f 


991  ♦ ***BBITE  (6,220) 

♦ 

♦ 

♦ 

+ All  80S  COMPLETES  ? (STEP  11) 


900  ♦ ♦♦♦♦♦♦♦CCNTI.IUZ 


CALCULATE  SUMBA81  RESULT  (STEP  12) 


♦ 00 

♦•*♦♦♦♦♦♦♦♦  27 

♦ ♦ xa»i,5 


(CONTINUED  ON 


pIgE  12) 


i 


A TOBPEEO  SIBOLAIIOH . 


bus  ?aoGa»s- 


? AGE  12 


♦ LB»1,IE0N 


♦ i 

♦ 

BB 

XH,  1 

♦ 

EB 

kb, 2 

♦ 

SB 

xa,3 

♦ 

SB 

xa.a 

♦ 

EB 

xa,5 

*SB(Xa,1  »ELOAT  (XCN (LB.XB)  ) 
»fB(XB,2  >DET (La , XB) 
*EB(KB.3  ♦ A 5S  j ASP  EC  ( LB  , X3)  ) 
*B8  (X3,4  ♦ 0ET3  (La,X3) 
c ‘"”=12, X 2 


= E B (KB,  5)  ► D EBE> 


28  ♦♦♦♦♦♦♦♦CCK1INUE 
♦ 

27  ♦♦♦♦♦♦♦ +CCKTINQ2 

I 

♦ 

♦ CO 

♦♦♦♦♦♦♦♦♦♦  920 

♦ ♦ XX*1,£_ 


C03POIING  BEAKS 

j 

| CEL  »5E(KX,1) 


* IE 
3 El  .12.  1. 


I DEL-1. 


53 (XX, 2 

aa  (xx, : 

EB  KK,  4 

ofl  XX,  5 


* S B (XX,  2 /DEL 
= EB  (XX,  3 /DEL 
»6fl  XX, 4 /DEL 
*BB  (XX, 5 /DEL 


920  ♦♦♦♦♦♦♦♦C05TIM0E 


(CONTINUED  ON  fiOE  13) 
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♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 


1 TOBFEDO  SiaOllIICN 


bum  Eaoafiia 


FACE  13 


♦—  ♦ 

♦ DO  ♦ 

25  ♦ 

♦ KK*1  ,5  ♦ 

j.Zl 

I 7AH<KK)  * 0. 

I 

♦ DO  * 

170  ♦ 

♦ 1*1, ISDN  ♦ 

"I 

. * 1 • . 

. • I?  * . 

• RCB  (I,**)  .fC.O  * 

* . . • 

* . . • I 

• . • 


I 170  | 


| VIR(KK)  * VI  B (KK)  * ( (DEI  ( I,  K K)  - 3H  (KK  , , 
I **2) 


) > I 


♦ ♦♦♦♦■♦ +CC  N I IN  OE 
25  ♦♦♦♦♦♦♦♦COSTIJIOE 


♦ DO 

♦♦♦♦♦♦♦♦♦♦  26 


♦ xa*i  ,5 


I CEl  *ai(Ka,1) 


. * li 
• DEI. IE. 2. 


. • 


1 DEL»2. 


| STD  (KB)  •SC3T(7an<Ka)  / (DEL-1.)  ) | 

' E3^3*1)  »?B(KB,  1)/FLOAT(I30S) 


i 

(CCSTIMOED  ON  fid  14) 


J 
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1 TORPEDO  SIEOLAIICN. 
♦ 

♦ 

♦ 

♦♦♦♦♦♦♦♦cchisiii 


bais  p£x>grab. 


PAGE  14 


EBIST  SOBHARY  (SHE  13) 

• ••MBITS  (6,  IS 7)  JR  ON 

FCEHA1  (1 X,//  ,61,' SOBHARY  0?  RESOLT  AFTER' ,3X, 14, 2X, 'RONS') 

• ••MBITS  <6, 19  9)  ( (RB  (I,J)  ,J=1 ,2)  ,STD  (I)  , (RB  (I,J)  ,J»3.5)  ,1*  1,5) 

PC E HAT  (101, ' tBCEA  EILIT Y OF  0 ETSCTIOK • , 3X . 


FCEBAT  (101,  ' r EC  c A cII.IT  Y OF  0 ETSCTIOK 3X . 

• AV'E RAGE', oX,' STD  DSVIATIOS'  ,7X,  ' AVERAGE1 ,7 X,  • AVERAGE*  , 
7X,' AVERAGE’, /,41X, 

■CET  RANGE' .BX,' DEI  RA MGE ', 7X ,' TARGET  AS? ECT • , 4X , ' D ET  EE 
3I,'R£I  EEAEXSG',/, 

IX, 'ONE  S0CCESSI72  D2T ECT ION 5X , F6. 4 , 3 J6X , F9 . 41 ,/, 


A&ISG ' , 


EGSBAI(1X,'HC  SET SCriOK  BADE  DO  RING  THIS  aOS') 


•£A2/iiCAT(IRiJN) 


•••WRITE  (6,  ISO)  C A 2 

?0EBAI(1X,/,  IX, 'AVERAGE  DEFLECTION  ANGLE  : • ,5X,F3. »,/) 
•••WRITE  (6, 2E4) 


• ••KBITS  (6. 236)  ( (DETS  (I , 1 ) , D2T  < I,  1)  , 1 :?SC  (I,  1)  , CLOSE  (1,1), 
CE1E  (1.2)  .DEljl.21  ,ASP»C  (1.3)  ,CLOsfl(I.i)  ,DSTS  (1,3)  , 

DEI  (1,3)  , AS  r EC  (1,35  , CiCS3  (1.3)  ) ,I«1,I3UN) 

FOR  HAT  (3  (IX,  56.  1 , IX,  76. 1,  IX,  ?6.  1 , 21,  F6. 1, 5 X)  ) 
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A TOE  P EDO  SIBOLATICB. 


SUSHOUTI'iE  PARSE!. 


PAGE 


I 


1 TOSEECO  SIHOLAIICN. 


SUBROUTINE  ? ABBET . 


S 0 E BO UTINE  EASHEI 
BZAEING  IB  DATA  AND  PABA3ETEHS 

COEBOB  ISEED2,  TTI3E,  TO,  TA,  T3AIE,  BANGE,  AL? A,  LA33D,  TADEC, 
EIAB,  B AC,  TAC,  CCOB,  DEVS? , SNG,  EN,  PH2,  SCCDHS,  ICOUBS, 

BIT,  HXS,  IK.  IDTI3E.  ITI3E,  XI,  IT,  ITAB,  ITAB , TDIST, 
aciol,  TUENTO,  INTVAi,  PHI,  B3AX(5,S),  TBANGZ , OISr.IPnlNI 

BE  1 1 LABEC,  HCOOBS,  BIT,  BXB , HOIST,  LA3BDG 

IBI2G23  SUNCUT 

TEDEC  - TECHNIC AI  DETECTION  BANGE  (STEP  A1) 

LEV  EL  C?  VARIATION:  315-150-1125-1500  3ETEBS 

TADEC  - TACTICAL  DETECTION  BANGE 

I 


f TEDEC*150. 

I IABEC*TEIZC 

I 


TUBE  - TBA  NSEISSICN  INISBVAL 


I 


I I II 3 Z“ 2. *1  ED EC/1500. 

I 


TCKN  - 10 Bp  SPEED  IN  KNOTS,  TO  - TOHP  SPEED  IN  3/SEC  (STEP  A2) 

LEVEL  CE  VABIATICN: 24-32-40  KNOTS 

J 

1 ICXB  *40.  | 

I TO  *TC  KN/2  | 


TAKN  - TARGEI  SPEED  IN  KNOTS,  TA  - TABGET  SPEED  IN  VSEC  (STEP  A3) 

LEVEL  C e V ASIATIC N:  12-13-24-30  KNOTS 

I 


IAKB  * 1 E . 

TA  *TAKN/2 


(CONTINUED  ON  PAGE  2) 


13U 


A TOE  FIDO  SlilULAIION 


SOSBO'JTINE  PARNET 


PAGE 


2 


I 

TACG  - TARGET  CCOBSE  IS  DEGREE  (STEP  A«) 

l 

I TACG  *270  | 

1 T AC  *TACG»RAD  ( 

I 

ALEAG  - SHEEP  ANGLE  IS  DEGREE,  ALFA  - SHEEP  ANGLE  IS  aADIANS 
(STEP  A 5) 

LEVEL  Cl  VARIITICS:  20-30-40  DEGREES 

[ 

i AiPAG*30.  I 

I ALFA  = ILF  AG* RAD ( 

I 

LAHBDG  - LOSE  HIETE  EACH  SIDE  OP  TORP  HEADING  (STEP  A6) 

IE7EL  Cl  VARIiTICN:  10-20-30  DEGREES 



| LA0BOG  *20.  I 

I LANED*LAGBCG*RAD  I 

I 

BELESG  - RELATIVE  3EARING  PROS  TARGET  TO  TORP  IS  DEGREE  (STEP  A7) 

LEVEL  CF  VASIJTICS:  0-30-60-90-120-180  DEGREES 

! 

I BELEfiG  *-60.  I 

I EEAR  * RELcEG *RAD  I 

I 

RANGE  - DISTANCE  EFTHEEN  TARGET  AND  TOBP  (STEP  A6) 

LEVEL  CF  VARIATION:  1500-3000-5000-7000  SETERS 

I 

I CANGE*  3000 . | 

I 

(CCSTINBFD  ON  pIgE  3) 
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A 


1 TC8PEE0  SIBOLATICS 


SOEBOOTISE  PAB3ET 


PAGE 


110 


I 

TBASGE  - 341  IOBP  fOS  IN  SETS3S  (STEP  49) 

J 

I I BADGE  -18CQ0.  I 

I 

IBAIEG  - TCBP  TOE  SEAT E ID  DEGSEE  PEE  SEC  (STEP  410) 

I 27 EL  OE  7 ASIATICS:  3-6-9-12-15-18-21  DEGSEE/SEC 

| 

I IBATSG  *18.  I 

] IBAT2=TBAT  EG*  BAD _ _ _ I 

I 

CAICOL ATE  WIDTH  CF  TACTICAL  SWEEP-LADE  (THEOBET ICAL) 

J 

| SB  DG  =TADEC*SID ( A1-PA+LA23D)  *2.  I 

~i 

CALCULUS  CCVESACE  BATIO  (TEEOBETICAL) 

I _ 

I CBAIIC  *1.-(TRAIE*TTiaS/(2. *14332))  | 


EEINT  CF  SITCATICN  AT  START  0?  SOD  (STEP  A12> 


• IF  < 

isaxsi.Ec.o 


. * 


I 90  | 


•**WBIIE  (6,110) 

PCE3AT(1X,//. IX, 'TACTICAL  3IIOATIOH  WHEN  FIBING',6X, 

' TCBP  EDO  EAEARETE5' 

21, 'BADGE  ATTACK  TARGET  T ABGET ' , 6 X. ' TEC . DET  TOBP' 
JX, 'SWEEP  XCSE  TOBD  SWEEP  C0VE3AGE' ./ . 

9 X , ' A DOLE  CCUBSE  SPEED', 7X, 'BADGE  SPafD',2X, 


•ANGLE  WICTE  SATE 
(CCNTIHOED  OS  pIgE  4) 


LADS' ,4X,  'BATIO') 
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k TORPEDO  SIMULATION.  S CESOUTINE  FIRING.  PAGE  1 

A TORPEDO  SIHOIAIION.  SUEBCUTINE  FIRING. 

SOEBOOTINE  FIRING 

CALCULATE  THE  IOBP  DEFLECTION  ANGLE,  BAIN  COURSE,  FIBING  CCDRSS 
BASED  CN  ESTIHAIE  OF  TARGET  DATA  (UNCERTAINTY) 


DIHENSION  0(E) 

CCBBON  ISEEE2  , TTI3E,  TO,  TA,  T3ATE.  RANGE,  ALIA, LA 3B  D ,T  A DEC , 
BEAR,  RAC,  TIC,  CCOR,  DEIS?,  3NG,  Pa , ?H2,  aCCURS,  TCOUBS, 
HIT,  aifl,  IK,  IDTI BE , ITIBE,  XI,  IT,  XTAR , XT AH,  TDIST, 

BDIST,  TOBNIC,  INTI AX,  PHI,  SBAX  (5,6)  , TRANGE,  DIST.IPSIMT 

CCEBON/DATA/Ci,OFLOS 

CCEHON/T AEG  ET/TAC3 G,TAB1 ,RNG8CD,DA1,COUR,CE (10)  ,SE(15)  ,JR0N, 
IE  I AG  , DA  2 

REAL  LABEL,  BCCUBS,  BIT,  BIB,  BDIST,  LABBDG 
INTIG  SB  RGNCCT 


PN  *-1 

1 

?P  * 1. 

K SPEED  *BCD ( 

;JRUN-D  , 15 ) ♦ 1 

1 

SCOURS  »IFIX 

(JBUN-1)/15.)  *1 

_J_ 

I 

CALL  GGOE  (ISIED2,2.U) 


CALCULATE  ESTIMATE  OF  TARGET  COUBSE  (STEP  31) 
1 

19  | TACB  »TAC*CE  (SCOURS)  I 

~~~ _l 

20  | CIECO-TACB-TAC  I 

I 

(CONTINUED  ON  pid  2) 
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1 T06EEEO  SIMULATION . SCBfiOOTINE  PiaiSG. 


FAG2  2 


. * * . 

. * If  • . 

* TAC3.GE.EB2  » 


( TAC3-TAC3-PB2 


. * • . 
. * IF 
* TACa.II.O. 
» . 

* . 

* . * 


r i TAca=PB2+TAca 


CALCULATE  ESI  OF  TA&G3T  SPEED  (STEP  32) 


j Till  *Ji*SS  (X3PSSD) 


CALCULATE  TOKFECC  CEPLECTICN  ANGLE  (STEP  33) 


| ASP  »EIA3-CTf CO 


. * * . 

. « IF  * . 

* asp.ii.-ebi  * 

* . . •- 

* . . * 

* . « 


| ASP«?H2fASP 


I ACS  *T  A S»  SIN  (AES  (AS  F)  ) 
CAA  »1 ACS/TC 


(CCNTINQ 


ED  CN  fLe  3) 


M 
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a to  bp  2 co  summon 


scebcutisb  Flams 


PAGE  3 


10 


. * 
* . 


• • * . 
IP 

CAA .Cl.  1. 


I 26  | 


| £ A = ArSIB  (CAA)  I 

I £A_  »SIGB  |EA,AS£)  I 

I 

CALCULATE  TO  5EE  DC  SAIN  rIEING  COUBSE  (STEP  34) 

I 


ECOUSS  »BBG*EHI*0A 

I 


- • * . 

» IP  » 

8C00ES.GT.PH2 


I 3COOaS=aCOOBS-?H2 


* IP 

accoss.u.o. 


I 9COOBS=P32+H COOES 


I — 


* IP  • 

ECOUBS.Gl .PH2 


* . • 


I 10  I 


(ccnimoPD  on  pIge  uj 


lliO 


1 


A TOBPEDO  SIBOL AIION . SCEBCOTINE  FIBISG.  PAGE  4 


CALCOLATE  TOEFEOC  PBES ENT.  FI3IJIG  COOBSE  (STEP  35) 


I 


i 

CIPALF 
ICOOBS  ■ 

'i0!1!*2- **!?»)  -ALPA 
'SCC  DBS* DIP  ALE 

1 

* 

. * 

. • i: 

ICOOBS 

1 

* . 

F * - 

.GI.EH2 

* 

* . 

* , 

. * 

. * 

r 

| TCOOBS*TCOOES-PK2 

t 

* 

. * 

• • i: 

ICOOBS 

* . 

F * . 

.LI .0. 

* 

* . 

. * 

T 

l TCoaas*PH2«-icoaas 

i 

1 

1 

* 

. * 

. * IJ 

DIPALE 

E * ' * . 
.GE.C. 

* 

* . 

• * 

T 

I EN=-PN 

i 

1 

CALCUIJT2  2S1IBAI2  OP 

TAfiGET  RANGE 

(STEP  36) 

• 

• * 

. * i: 

0(2)  .( 

i 

* . 

p * . 

51.0.5 

* 

* . 

• * 

I 

1 PP—PP 

1 

i 

! 

| B1IG0IP  » (1.-0  C2) ) *BAHG!*0. 15  | 

(CCNTI80E0  ON  P, 

Ice  5) 

llll 


""" 


A TOB PiCO  SIHOL ATICN.  S OEECOTIN2  PIBING. 


JUGBOD-BASGE+SIGN  (BNGDIP.PP) 


IACBG*TAC3,/BAD 
TIBI  *1  ACS  * 2 
EA1  -CA/BAC 
CA2  »Ci2*CA1 


I 


PAG2 


122 


124 


125 

126 


PHIHT  COT  OP  PIBING  DATA 

. * IP  •■  . 

* IPBISI.EC.1  * 

• . . 

* . . * T | 25  | 

♦ * ” ————— — 


***KBITE<6,  122) TACBG, IAal.BNGBOD 

POBMATMX,  'EST  OP  TAP.G2T  DATA  ?OH  PIBING',/, 
4X,' CODES* ',SX,' SPEED' ,61 , ' B AHGE' ,/, 

1X,J<?8.  1,32)) 

•••KBITS  (6 , 124)  C A 1 

PCBHAT(1X,'TCEE  DEELECTION  ANGIE  IS  ',P6.2) 

| COOa  • BCC  G BS/B AD  I 

I 

•••WHII2  (6,  125)  COOH 

pcaaAnu.'TcapEoo  aiis  coubse' ,ax,F6. 2) 

CC  N UNDE 

I 


26  •••KBITS  (6 , 2C) 

30  PC  E NAT  (/  , IX, ' NCT  P E AS  ISLE  TO  FIB5  COE  TO  LOW  T0BP2D0  SPEED') 

| IPLAG»  1 | 

*~l 

25  BEIOBN 

ESC 


11:2 


a 


A TOFEEDO  S ISOLATION . SCESOOTINE  FOSIS. 


1 TCSEECO  5 ISOLATION.  SOEHOOTINE  FOSIS. 

30ESO0TI&E  ECSIS 

IS  CALCULATING  NEK  POSITIONS  OF  T ASSET  AND  T09PEE0  IN 
EACH  TIBI  STEF 


COB  BON  ISEEC2,  TUBE,  TO,  TA , TBATE,  3ANGE,  ALFA,  LAK2E,  TiDEC, 
SEAS,  .U[,  TAC,  CCCS,  C2VSP,  SNG,  tS,  PH2,  3CC0SS,  TCOUaS, 

BIT,  BIB,  IK.  I DTI  BE,  ITIBii,  XT,  IT,  ITAS,  TTAE,  TDIST. 

BCISI,  TufiNIC , INTVAi,  PHI,  3«AI<5,o),  T3ANGE,  5isT,I?bINT 

BEAL  LABES,  SCCU3S,  BIT,  BIB,  BDIST,  LABBDS 

INT2G23  SONCOT 


TIBE  COUNT  (STEP  Cl) 


| IK  »IK+1 
I IIIBE^IIIBE+IDTISE 


CALC3IAI2  TOTAL  ICHP  30N  AND  TABG2T  BUN  [STEP  C2) 

I 


| aiT  »elT*TIIST 
I exa  =kxb*sdist 


CALCULATE  NEK  POSITIONS  (STEP  C3) 

I 


* XT  ♦ SI  N (TCOURS)  «TDIST 
«TI*CCS JTCOUKS  *TDIST 
SIN  (TAC)  *f 


HAS  3XIAS*SIN  (TAC)  «SDIST 
HAS  3YTA3+CCS (TAC) «BDI5T 


CALCOLATE  NEK  TOSE  C008SE  (STEP  C4) 


( IXC0U3  *TCC  D8S*SIO-N(TU2STO,EN) 
I 3 A E S (BCOUSS-IXCOUa) 

(CONTINUED  ON  fIgE  2) 


A TOBPECO  SIMULATION 


SCEBOUTIITE  FOSI5 


PAGE  2 


A TOBFEDO  SlitJLAIICN 


S CE BOUTINS  DETECT 


PAGE  1 


A TOBPEDO  SIECLAIICS.  SOEHOOTI ME  DETECT. 

SDBBOUIINE  DETECT 

TO  CHECK  I?  T1BGZI  IS  DETECTED  AND  STOEE  DETECTION  DATA 


CCHHON  ISEED2.  TTI3L,  TO,  TA,  TEATS,  HANGS,  AL/A,  LA9BE,  TADEC, 
SIAB,  BAD , TAC.  CCC5,  DEVS?,  3NG , BN,  ?H2,  SCCUBS,  TCCOSS, 
an,  sxa,  ig.  idtise.,  itiss,  xt,  it,  xtab.  stab,  tdisi, 

SDIST,  TOENTC,  INTVAL,  phi,  H£AX  (5,0)  , TSAMGE,  DIST.IPEINT 

CCEflCN/DATA,/H,OPLCE 

BEAD  LAflED,  SCCUBS,  HIT,  9X9,  HOIST,  LAH3DG 
INTSGEB  fiONOCT 
DIB  ENSIO  N D 8 (3,2) 

DOUBLE  PBECISICN  PCWB AX, 3 , R3X , S3 1 1 , 38X2, X 1 ,X2 ,X3 , XX 1 , 
V1,V2,V,U,fIE ACT, X Cl ,X0i,X03 ,XX2,?0MSB,3£13,AiAa 


SETTING  0?  TABG IT  DI3E3SI0N,  A 


IA3GET  LENGTH, 


B - TAEGET  HIITB,C  - TABGST  DEPTH.  (STEP  D1) 
I 


1 A 

» 100. 

t 

E 

- 15. 

i 

1 c 

»4. 

I 


CHECK  II  T 9 A N53ISSIGN  (STEF  D2) 


I 

. • • . 

* II  • 

IK.LT.INTVA1 


I 20  I 


I IK  »0 
(CONTINUED  ON  I 


la  £ 2) 


1 TOSPZDO  SIfULAIICN 


SOESCOTINS  DETECT 


PAGE  2 


CALCULATE  iAKGE  TC  TABGST  (STEP  23) 

I 


CIPX  -XTAE-IT 
DIPT  *XIAH-IT 

CIST  »SCST  | (DIPX»*2) ♦ (DIPT*^) ) 


TEST  IT  TAflGZI  IS  SITHIN  POSSIBLE  3ETECTIO N 3ANGE 


. * I?  * . 

■DISI.GI.  (TA2 ZC+A/2  . ) 


I 20  | 


DETECTION  TH3ESBCLE  (STEP  2U) 

J 

ECNHAX  *1./  (B*TA32C)  **<* 

U •‘SCALE  If  HI/ 2.  | 

CONST*  <A**2)  • (3»*2)  • <C*< 
rOUNAX  *CC  N ST*PO  NNAX  *(J 

I 


2) 


CALCULATE  TLNE  TC  TAflG ET  (STEP  25) 

I 

TIS2L1  *0131/1500, 
XTAS1»Xli3*SIN  (I  AC)  • (7I«DL1*TA) 
YTAB 1«ITAS»CCS  (TAC)  * (ZiaDL1»TA) 

f 


CENTS  - BEASING  CP  C2NTE8  CP  LOBE  (STEP  06) 

j 

I CD  *-CA»CPLOB 
| CENTB»TCCOSS*CD 

I 

(CONTINUED  ON  fid  3) 


lu6 


1 TOBPZEO  SleULATICN.  SDEEOOTINE  CETECT. 


PAG2  3 


. * * . 

. * IF  * • 

• C2NTE .II . G.  « 

• . . • 

* . . * 

* . * 


| CENTS* PH2+CENT5 


. * * . 

. » IF  * . 

* CENTS. GI.EE2  * _ 

* . . * T 1 CENT B*C2NTB-PH2 

* . • 


CIIC0IJT2  SEAFISGS  TO  TA3GIT  (STEP  07) 


•EEAFIN (XTA31 ,YTA51,IT,YT) 
= 2Ii3  1*:=IN  (TAil)  »A/2. 
s 1 1 A H 1 +COS  TACi  *A/2. 

•SCSI HXIAB2-XT) **2» (ITAH2 

• EElllN)j(TA32.YTA32«XT,YT) 
•XIAB  1-EIN  (TAO  *A/2. 

•Yiiai-cos  (TAc5  *a/2» 

•seal  I UIAB3-XX) **2» ( Y1AB3 
-»T)  *»2) 

•EEASIN  (TTAB3,ITAa3,XT,YT) 
•BE  1— CENTS 


. * • . 

• IF  • 

BEL1.GT.FBI 


B2L1=BEI1-PH2 


(CCNTINDEO  ON  PAGE  4) 


Hi  7 


A TOBEEDO  SieULAIICH 


SOEIOOTIH2  CEXECI 


PAGE 


4 


. * * . 

. * II  * . 

• BIL1.L1.~E3I  • 

• . ^ * . * T | BEL 1»PH2+  BEL  1 


1 BELO  =Sf-CESXB 

1 

1 

. • * . 

. * IE  * . 

• aSL0.C-l.E3I 

* 

* * . . 

* . * 

;■ 

r | BELO2  BELO- PB2 

1 

. * * . 

. * if  * . 

* 3EL0.L1.-EHI 

* 

* • , * 

* _ * 

1 

X | BELO2 P 32+ BELO 

i 

| 

1 EEL2  »a£2-C£NTB  1 

1 

1 

. » * . 

. * IF  * . 

* BSL2.GI. FBI 

* . 

* . . * 

* . * 

* 

* 

X | BEL2»B£L2-EH2 

(CCHXIHOZD  OH  ?1gE  5) 

A T03PED0  SINOLATICN 


S OE  f CUT  I N £ DETECT 


?4GE 


5 


I 

. * * . 

• * IP  * 

* 3IL2.LI.-E8I 

• . 

* . . • 
* . * 


| HEL0  = PH2*  3EL2 


CALCULATE  TBANSNISSION  GAIN  EACTOR  (STEP  08) 

^ 

33X1  =3111 
SEX  = 5110 
EBX2  =3112 
ALAN  = L AN  £ E 


CONFUTE  SEPASAT1  GAIN  2 ACT  CHS  (STEP  09) 

I 


1 11 

=Xf ACT 

(33X1 

1 X2 

=X1 ACT 

(aaxx 

! 13 

=X1 ACT 

lSBXf 

XXI 

= (inx 

J + X3) 

CALCULATE  TA3GET  ASPECT  ANE  TA3GET  SON A3  C30S5- 

A3EA  (TiBGET  STS1SGT9)  DUE  TO  ASPECT  (STEP  010) 

I 


| A NTH3=  5 1 *E  El 


. • 

• II  » 

ANia3.CT.EB2 


* . 


* 


| AMTB3= ANT3B-EH2 

I 

/ 


| SELA  * A NTS  1 — I AC 
N e1g£ 


I 


(CCSTISOED  ON 


6) 


1L9 


A TOEPEBO  SIHULAIICH.  SCEEOUTIME  DETECT.  PAGE  6 


. * If  * . 

• aSLA.GI.EHI  * _ 

* . . * f | B3LA-PH2-8ELA 

* . * 

J I 

i 

. « * . 

• » If  « . 

• 3ELA.L1.-Eai  J 

• . . « I | B£LA=PH2+ BELA 

*i‘  “V 


5 EL  £ 
VI 

V 2 

V 

=*AES  (BELA) 

= (A*«2)  * (bcCS(5EL3) 
= (E**2i  * (D5IM.(EEL£) 
» (V1  + V*)  »»2 

**  2) 
•*2i 

1 

! 

i 

1 

COaiCTI  SCAXISG  c ACTOH  DUE  TO 

1 

ASPECT 

(STEP  Oil) 

1 U 

1 f X?A 

= EC  A LE (BEL3) 

CT  =0/V 

1 

1 

l 

CALCULATE 

1 

BETUBNTIBE  ?CE  ECHO 

1 

(STEP 

312) 

1 TIHDL2  =*2.*TiaDL1 

1 

CALCULATE 

1 

BEL  EEASIHG  ?CE  SETUaHING 

1 

ECHO  (STEP  D 13) 

1 TOBMST  =THATE«TiaDL2 
| TIC  *TCO OSS ♦SIGH  (TUSXST 
1 1 I3C  =AES  (HCCUSS-T2C) 

,?X) 

1 

I 

I 

(CCXIIMOEB  CX  fief  7) 


150 


A IOEFEDO  AiaOLAIION 


SOEBOOTISE  DETECT 


PAGE  7 


« IP  * 

TXDC.CT. FBI 


_ m 


| TXDC=PH2-TXDC 


* If  * 

TX2C.IE.AIf A 


18  | 


I ALP  -T  XL  C *•  AL*  A 

I BELT0=5EL0-EIGN ( ( TUB NST- ALP)  ,?N) 
| + SIGS  (ALE,  EN) 


I 19  | 


18 


t EEIT0=PEI0-5IG8(TUENSI,PM) 


19 


3EIDIP  = 5 EXT  C- 3E10 
BELT  1 = SEL 1 *62X21 r 
3ELT2=5EL2+3ELDIf 


1 


- » • . 

. * If  * 
« PEL  10. €1. PEI 
* . 


* . * 


| a£LT0=PH2-EElI0 


I 


I 


(CONTINUED  ON  EAGE  3) 


A TOE  PECO  SIKULATIOH . SCEBOUTIHg  DETECT. 


PAGE  8 


* * . 

* I?  * 

8ELI0 . II . -EHI 


I BELT0-PH2+BELT0 


. * rr  • . 

« S21T1.II.-EBI  * 

* . . *- 
* . . * 


T | aELT1=PS2+BELT1 


. * IE  • 

» BELI2.GI.EHI 

* . 

* . . « 


* 


: | BELT2-PH2-BELT2 


. * * . 

. • IE  * 

* BELI2.II.-Eai 

* . 

* . . « 
* . * 


| BELT2»PH2*3ELI2 


(CCHIINDED  CD  p la  £ 9) 


. • IE  * 

» aELI1.Gl.EHI 

* . 

» . . * 
* . * 


T | BELT1=PH2-BELT1 


A TOE? EDO  SIMULATION 


SOEBOOTINE  DETECT 


PAGE  9 


TEST  ECS  DOPPLEfi 


EBQDBF  »50.*2.*(T0*CCS(ED)-1.) 

/ 150C. 

EBQDIF  »50.*2.«(TO*CCS(DD) *1.) 

/150G. 

EELC  *5£LE 

EBQS8-50. *2 . * (TO *COS (DD) ♦TA*COS (8ELC)  ) , 
| / 1500. 


* (EEQSa.il. 7BCDIE)  .AND.  (F2Q5H.GT. FHCDEF) 


CALCULATE  2ANCE  AND  3EA2ING  TO  CLOSEST  PAST  OF  TABGET  (STEP  D1U) 


1 3 

i 

1 

•DIST  1 

:a 

i 

2 

•2EIT0+0D 

ca 

2 

1 

=C IST1  1 

C 3 

2 

2 

«aELX1*0E 

e s 

3 

1 

•DISI2  I 

EE 

3 

2 

=»aELT2*0E  | 

cai 

N 

= 

:b(i.i)  I 

ao 

3 

. * IE  * . 

* CaiN.GI.E3  (2,1)  * 


I 10  ( 


» EaiS.Gl.E3  (3, 1)  * 

* . . »- 


111  I 


(CONTINUED  ON  PAGE  10) 
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A T08P2D0  SIBULAIICH 


SCESOUTIH2  DETECT 


PAGE  10 


10 


11 

12 


I 


| CBIX  =Ll  (2,1) 
I BE  3i 


I 9 | 


I 

i 


ESIX  = CE  (2,1) 
BE  =3 


I 


| ESEL  =CE  (BE, 2) 


* 

* . 


* IE  * 

3SEL.GI.Pai 

« . . * 


* IE  « 
BS  EL  . LI E HI 

* . . < 


BS2L=PB2-3SEX 


| 3 8 EL3  PH2  ♦ 3 B EL 


CALCULATE  SECEITISG  GAIH  EACTOB  (STEP  E15) 
I 


8 B X 1 

»8I 

!LT  1 

sax 

=>s 

-LIC 

S 3 X 2 

■B 

ILI2 

X 0 1 

»x. 

ACI  IBEX  1 - ALAS) 

X 02 

3 X 

ACI  (BSX.ALAS) 

103 

«x: 

ACI  (83X2, ALAH) 

XX2 

* ( 

C 1«I02*xi31  /3.  1 

(CCSIISOEE  OH  EiGE  11) 


l$li 


A 


A IOSEEDO  SIBUUIICH 


5 CE BOUTINS  DETECT 


PISE  11 


I 

CALCULATE  FSACTICN  OP  POMEB  IN  TO  3ECIEVZB  (STEP  016) 

I 

| E0HEB*CCNST*I11*XX2»PI?ACT/  (DIST  | 

I 


TEST  FOB  DETECTION  THHESKOLE  (STEP  D17) 


. » • . 

. * I?  * . 

* PONES. IT. PC » SAX  • 


COHIOII  EEABIBG  !ATE  (STEP  D1 3) 

I 

( ASP  =8S*PHI-TAC 


« TV 

ASP. ST. PEI 


ASP-ASP-PH2 


. • 

» IF 
ASP.LT.-EHI 


• . * 


( ASP’ PH2*  AS  ? 


| AP  = ASP*  BELTO 


(CCNTINDED  ON  PAGE  12) 


I 
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1 


i TOBPECO  SIKULATICN.  SCESOUTINE  DETECT. 

I4CS«T**s|s  (iES (ASP) ) 

T'iocs  «TC*S I S ( ABS  B2LI0)) 

I cBAIE-  (IACS*SIGN  TOCS. A?)  ) /DIST 


CHICK  BEARING  Hill  AGAINST  T0RN3ATE  (STEP  D19) 


PAGE  12 


. • * . 

. « IE  * 

* EHATE.GS.I3ATS 

• . 


I 15  | 


CHECK  IOBPEDC  SPEED  AD  V AST  AG  2 (STEP  D20) 


1AIS  *TA*CCS(A23  ASP)) 
IOIS  »TC»CCS  (A3S  JaEIID) 


. * * . 

. * IE  * . 

•AES  (ASP)  .GI.  PHI/2.  • 

* . . *- 
• . . * 


| TALS»- TA1S 


. » • . 

. • IE  • . 

•AES  (aELTC)  .GI. PHI/2. 


| TOLS»-TOLS 


(CCNIINOID  ON  PAGE  13) 
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1 TOSPEDO  SI E OL ATION 


SCcSOUTINE  DETECT 


PAGE  14 


xcosi.Gr.s 


ail 


A T0SP2C0  SIBOLATICN 


SCEBOUTISB  DETECT 


. * X?  « . 

* EBAX  (2, 1)  . SE.O.  • 

* . . *• 


* . 

* , 

. * iT 

, * 

1 

B3AX  (2,1  »DIST 

EaAX(2,2  * (BELTO-KJD) /8  AD 

. 

. * 

. * 11 
j aaiX(2,2) 

* . 

* , 

* . 

* 

'.Gi.iaol  * 

• * T 

• 

• 

1 bbax(2,: 

| EaAX(2,i 
I SaAX(2,I 
saAX(2,e 

i)  =»DBIN 
>)  *ESEX,/BAD 
:)  *SELA/SA0 
:)  sSLi/3AD 

1 2C  ( 

.J 

. » ii 

* SBAX(3,1 

* . 

* . 

* . 

* . 

1)  .NE.O.  * 

. * T 

* 

i 

S8AX{3,1>  »DIST 

EB AX (3,2  * (H£iTO*OD)/HAD 

(CONTINUED  ON  fi 

G E 16) 
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A TOEEEDO  SIEOLATICS 


SCEEOOTIHE  DETECT. 


PAGE  16 


. * * . 

• • If  » • 

» sail  (3,  2)  .Gl.  180.  * 

* • . • *“ 
• . . » 

* . * 


| BBAX  (3.2)  *8BAX  (3,2) -360. 


B B A X i 3 , 2 
SBAX  (3.4 


=cais 

=£ BEL/ SAD 


EMAX(3,5J  ’BEIA/HAD 
E 3AX  (3 ,6)  =El£/3AC 


I 20  | 


. * * . 

. * IE  ■ . 

* EBAX  <4 , 1)  . SE.O.  * 


I 20  I 


| BBAX(4,1)  ’CIST 
| 3BAX  (4,2)  * (RSLI0*DO)/aAD 


. • IE  * . 

* SB  AX  (4,2)  .GT.  180. 


| HOAX  (4,2)  «aSAI  (4,2) -360. 


EBAX  (4,3  ’CflIH 
E B AX  (4 , 4 ’EBEL/EAD 


EBAX  4,5)  ’8ELA/3AD 
EBAX  4,6)  *ELE/'aAD 


(CCHTIHUED  OS  PAGE  17) 


A TOBf EDO  SIEOIATICN.  SCEBC0Tl;i2  DETECT. 


PAGE  17 


I 20  | 


• £ a AX  (5,1)  . BE. 0.  * 


I 20  | 


BHAX(S,1  =OIST 
BHAX(5,2  * 18ELTO»OD)/HAD 


. * * . 

. * If  * . 

* EH  AX  (5, i). 01.180. 


I BHAX  (5,2)  =BHAX  (5,2) -360. 


I 3H  AX (5,3  *CHIN 
BHAX  5,4  -cS2L/ BAD 
I BHAX (5,5  -BBLA/BAO 
I BHAX (5,6  “BIc/BAO 


IF  SO  CETECTICN,  BET  DETECTION  3IAT3S 


ICO  NT*  0 
JCONT’O 


A TOBPEDO  SIBOLAIION 


POSCTIOH  BEABIN 


P AGS 


1 


16 


I 

A TOBPEDO  SiaOLAIICH.  FUHCTIOH  2EABIH. 

FOHCTIOH  etioIS  (i,a,C,D) 

10  CALC OL ATE  EEABIHG  EBCa  IOBPSOO  10  TABGEX 

I 


I CIEX 

*A-C 

| 

| LIE  X 

»a-D 

I EH2 

= 2. *3.  141  = 92654 

| 

1 i AO 

*EE2/ joO. 

i 

"1 
. ' • . 

. • 

IE  * . 

* DIET. HE. C.  * 

* . . * 

« . . * T | 16  | 

’ i ‘ 

i 


| EE  »90.*EAC | 

. * » . 

. * If  * . 

» DIE X .LI. C.  • _ 

* . . « X | H3=«B+  (180 . »SAD) 


I 17  | 


_i 

I BS =AIAH2  (D I?X,DIET) 

(CCHTISOEO  OH  eLe  2) 


162 


I 


A TOE  BEDO  5ISULAII0N . EOHCTIOS  3EABIS. 


PAGE  2 


17 


ESC 


| E £ a aiti  =aE 

I 

azicas 


i 


163 
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A tO  E p £ CO  SISOtATICH.  FUNCTION  XFACT. 
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